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1990 Workshop Summary 


The 2nd NASA CFD Validation Workshop was held at Lewis Research Center on 
July 10-12, 1990. 

The purpose of the workshop was to review NASA's progress in CFD 
validation since the first workshop (held at Ames in 1987) and to affirm 
the future direction of the NASA CFD validation program. , 

Forty-six people participated in the workshop. Of these, 25 were from 
NASA, 2 were from other government agencies, 13 were from industry, and 6 
were from universities. A list of attendees is included in this report. 

The first session,, held on the morning of the first day, consisted of 
overviews of CFD validation research at each of the three OAET research 
centers and at Marshall Space Flight Center. The second session 
(afternoon of the first day) consisted of in-depth technical presentations 
of the best examples of CFD validation work at each center (including 
Marshall). 

On the second day the workshop divided into three working groups to 
discuss CFD validation progress and needs in the subsonic, high-speed, and 
hypersonic speed ranges. The emphasis of the working groups was on 
propulsion. The subsonic, high-speed, and hypersonic working groups were 
led by Gordon Pickett (Pratt and Whitney), Joe Marvin (Ames), and Lou 
Povinell i (Lewis), respectively. At the end of the second day, each group 
leader reported to the workshop on his group's findings and 
recommendations. The rosters of the working groups, and copies of the 
charts used in reporting out are included in this report. 

It was apparent from the presentations that NASA's CFD Validation program 
has significantly changed the way experimental fluid dynamics research at 
NASA is performed - even research that is not funded out of the CFD 
Validation Program. - Compared to before the validation program began, 
researchers today are finding it easier to obtain funding for 
instrumentation and tunnel time required for code validation. As a 
result, more experiments are producing data that is useful for code 
validation. The synergism between experiment and computation is producing 
better experiments and better computational methods. | 

CFD^Val idation Coordinating Board Recommendations | 

On the third and final day the NASA participants met to discuss the future 
of the NASA CFD validation program in light of what we had learned in the 
first two days. The Board, with input from other NASA participants, made 
the following general recommendations: 

1. NASA has made significant progress in implementing the recommendations 
of the Ad Hoc Committee on CFD validation of the Aeronautics Advisory 
Committee and should continue to abide by them.. 
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2. NASA has made significant progress in implementing the recommendations 
of the first NASA CFD Validation Workshop and should continue to follow 
them. 

3 Changes in the work breakdown structure for FY91 and beyond must not 
adversely affect funding for CFD validation. In particular, since with 
the new WBS CFD validation will no longer have a funding stream 
separate from that for CFD Methods, NASA Headquarters and the research 
centers must ensure that the CFD validation programs is adequately 
funded. 

4. NASA must maintain a highly visible CFD validation program to enhance 
advocacy for obtaining funding from other sources. 

5. NASA should encourage applied programs to incorporate into their test 
plans methodology that allows for some level of code validation. 

6. NASA should advocate to program offices that they provide advanced 
instrumentation to quantify flow environments. 

7. NASA should encourage the continued participation of Marshall Space 
Flight Center in the validation program. The Board commends their „ r 
synergistic role in working with the research centers. 

8 The Coordinating Board should organize a third NASA CFD Validation 
workshop to be held at Langley Research Center in approximately three 
years. 

The Board identified the following important propulsion validation 
needs that are not being adequately addressed: 

1. supersonic mixing and combustion (2<M<6) 

a. lower speed: instrumentation is lacking 

b. higher speed : faci 1 i ties are lacking 

2. subsonic combustion 

3. transonic afterbody propulsion integration 

4. secondary flow systems in propulsion systems 

5. supersonic/subsonic jet mixing for HSCT 

6. swept-wing aircraft icing 

The Board recognized that, since not all cognizant people were in 
attendance, this list may not be comprehensive. 
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April 23, 1990 


Dear: 

I would like to invite you or your representative to participate in the 
Second NASA CFD Validation Workshop, to be hosted by the Lewis Research 
Center, July 10-11, 1990. 

In 1986, NASA's Aeronautics Advisory Committee (AAC) formed an Ad Hoc 
Committee to review the CFD validation activities at the three NASA 
Research Centers. The Ad Hoc Committee was charged to: 1) classify CFD 

research to identify and characterize verification needs and, 2) assess 
ongoing and planned CFD verification experiments. The ad hoc committee 
conducted reviews at each of the three NASA research centers and reported 
their recommendations back to the AAC. In response to those recommenda- 
tions, NASA conducted the First NASA CFD Validation Workshop, held in the 
summer of 1987 at the Ames Research Center, as part of the process to 
insure that the focus and future direction of the validation program was 
consistent with the anticipated needs of the aerospace community. 

The time has now come for a follow-up second workshop. The purpose of this 
workshop will be to summarize progress in NASA CFD validation experiments, 
instrumentation, and facilities since the first workshop, and to affirm the 
future direction of NASA's CFD Validation Program. 

It is planned that the first day of the workshop will consist of 
presentations summarizing the progress in CFD validation at the three 
research centers and at the Marshall Space Flight Center. The second day 
will consist of small working group meetings, focusing on specific topics. 
Since we will be meeting this year at the Lewis Research Center, the 
emphasis of the meeting will be on propulsion, focusing on the three 
principle speed ranges - hypersonic, high-speed, and subsonic. 

Your active participation in this workshop will be of great value to NASA 
in insuring that the CFD Validation program continues to make progress in 
the areas of most critical need. 

Please return the enclosed registration form as soon as possible. A 
detailed program, information on hotels and direction to NASA Lewis will be 
sent to registrants at a later date. 

If you have any questions, please contact the Workshop Chairman, Dr. 

Raymond E. Gaugler, at (216) 433-5882. 

Sincerely, 


Neal T. Saunders 
Director of Aeronautics 
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Conclusions and recommendations of the 
Advisory Committee, May 1987 


COMPUTATIONAL FLUID DYNAMICS VALIDATION 


Prepared for 

NASA'S AERONAUTICS ADVISORY COMMITTEE 

by . 

THE AD HOC COMMITTEE ON CFD VALIDATION 
OF THE AERONAUTICS ADVISORY COMMITTEE 


May 1, 1987 
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CONCLUSIONS 


Many excellent experimental programs are bei ^ ejected at 
«. "ci rpsparc h centers. The programs are generally attaching 
Srt inent problems . However, only a few of the projects reviewed 
were^ considered by the committee to be C™-velidation project^ 
Further adequate coordination among NASA centers was nor 
evidenced in the reviews. 

CFD validation is a relatively new concept that requires 
clearly defined and disciplined coupling between CFD and experi- 
ments/ The committee has proposed the following 
CFD validation which emphasizes the computational and expen 
SSJl ^discipline required end distinguishes velideticn from 
other experimentation related to CFD: 

rv n code vaUdation . - Detailed surface- and flow-field 
comparisons with experimental data to verify the code s 
ability to .irately model the critical physics of the 
flow. Validation can occur only when the accuracy and 
limitations of the experimental data are kn°^. 
thoroughly understood and when the accuracy and limitations 
of the code's numerical algorithms, grid-density effects, 
and physical basis are equally known and understood-* over a 
range of specified parameters. 

The committee recognizes that four categories of experimen- 
tation are required for developing CFD capability: 

A. Experiments designed to understand flow ph ysics. 

B. Exper im ents designed to develop physical models for CFD 
codes . 

C. Experiments designed to calibrate CFD codes. 

D. Experiments designed to validate CFD codes. 

All four categories of tests are important and are necessary 
to build a mature CFD capability. Validation tests shouldbe 
only a part of total test focus of NASA and should be formulated 
to provide specific data for validating CFD codes. 

CFD validation is severely hampered in some areas by the 
lack of critical measurements under realistic conditions- 
especially at high Mach numbers. Measurements must J be taken with 
adequate accuracy and resolution, and with redundant instnments, 
in order to evaluate CFD's capability to predict details and 
trends as well as to explore boundaries of application for 
specific codes. Specific instrumentation and facilities need to 
be developed when the state of the art is inadequate. 

Equally important for validation is the mapping of the CFD 
code's sensitivity to numerical algorithms, grid density, and 
physical models. These effects must be known to the same degree 
as Y the experimental accuracy and resolution in order to under- 
stand the applicability of the code over a range of flow para- 
meters . 
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RECOMMENDATIONS 


It is recommended that 

l NASA classify its CFD related experiments into the above- 
defined categories A, B, C, and D. The classification 
provides a framework for future evaluation and allocation of 
resources . 

2. NASA adopt a sharply defined program plan for CFD valida- 
tion, consistent with the definitions proposed by the 
Committee in this report. In particular, 

o Projects falling under this program should be closely 
coordinated and adequately funded. 


o The limitations of the facilities, instrumentation, 
proposed models, and codes should be clearly specified at 
the onset of the projects, as should the critical data 
set required and the range of validity of the modeling 
assumptions. 

o The questions of such phenomenon as unsteadiness, three- 
dimensionality, transition, turbulence models, and the 
applicability of inviscid/interacting-viscous/Navier- 
Stokes models should influence the choices of experi- 
ments, code development, and application. 

o Duplication of experiments and code development should 
involve several centers . '■ 

o A validated code should become a NASA or general code, 
not merely the code of a given center • 

3. NASA pay careful attention to the instrumentation and faci- 
lities required to provide data of adequate accuracy and m 
the form and detail required for CFD validation. Close 
coordination must be maintained between centers. 

4 NASA establish a Coordinating Board for CFD validation, 
consisting of key project leaders from each of the three 
NASA centers and chaired by NASA Headquarters . This board 
should have responsibility for defining, coordinating, and 
focusing the CFD validation effort. 
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WORKSHOP SUMMARY 


The first NASA CFD Validation workshop was held at the Ames Research Center 
on July 14-16, 1987 with 108 persons in attendance. After being welcomed 
by the Director of the Office of Aerophysics at the Ames Research Center, 
NASA Headquarters management presented introductory comments and expecta- 
tions for the meeting. Recommendations from a recent review of CFD valida- 
tion activities at three NASA Research Centers, conducted by the Ad Hoc 
Committee on CFD Validation of the Aeronautics Advisory Committee (AAC), 
were presented. Definitions developed by the Committee for code verifi- 
cation, code calibration and corresponding experiments were also introduced 
and discussed. NASA will use these -definitions in conjunction with the CFD 
Validation program as a tool to manage and maintain a balanced program. 
For completeness, a copy of the AAC Ad Hoc Committee report is included 
herein. 

The Ames, Langley, and Lewis Research Centers presented overviews of their 
CFD Validation activities with brief descriptions of their key experiments. 
These presentations covered the entire speed range and included external 
and internal aerodynamic flows, combustion, aerothermal loads, and aero- 
elasticity. Following the NASA presentations, other government., instal- 
lations, industry, and universities presented overviews of their CFD 
Validation activities and needs. The presentations and speakers are listed 
in the Workshop Program which is included herein along with copies of the 
presentations. 

After these formal presentations, the workshop attendees were divided into 
six working groups. The working groups focused on the fallowing key areas: 
Low Speed Aerodynamic Flows; High Speed Aerodynamic Flows- Commercial; High 
Speed Aerodynamic Flows-Military; Internal Flows; Hypersonics-Chemically 
Reacting External/ Internal Flows; and Hypersonics-External Aerothermo- 
dynamics. Each working group was given a list of five issues to address in 
their discussions. The issues were: (1) identify key efforts (numerical 

and experimental) required to meet immediate modeling and validation needs; 
(2) identify near-term and far- term critical problem areas that require new 
or additional modeling and validation ‘ activities; (3) identify computa- 
tional, experimental, facility and instrumentation, or other capabilities 
required to investigate these critical problem areas; (4) identify key 
modeling and validation projects that ' are potential cooperative/joint 
ventures; and (5) define, identify and/or propose standardized test cases 
for modeling and validation. 

The working groups met for about 6 hours to discuss these issues and 
prepare summaries. 

The formal workshop resumed with the various spokespersons presenting the 
results from their group discussions. These presentations have also been 
included herein. Several common recommendations were made. These, not in 
priority order, included: (1) provide closer cooperation between CFD 
developers and experimentalists at the outset of verification projects with 
a lasting commitment from both to see the projects through to completion; 
(2) provide detailed measurements of the flow field and boundary conditions 
in addition to model surface measurements and integral quantities; (3) pro- 
vide improved or new non-intrusive measurement capabilities, especially for 


hypersonic or reacting flow conditions; (4) provide redundancy in both 
measurements and experiments whenever practical so as to clarify data accu- 
racy and credibility; (5) provide dedicated large facilities for validation 
research and increase flight-based research activities; and (6) provide 
standardized test cases with accessible electronic data bases. The recom- 
mendation for standardized test cases received a strong consensus during 
the workshop. NASA will review the recommendations from the groups and 
determine a course of action to advocate their implementation. 

The workshop concluded with the general feeling that the meeting was very 
successful. Attendees from the various research centers, universities and 
industry expressed a renewed understanding of each other's viewpoints on 
CFD validation. Experts in computational and experimental fluid dynamics 
were able to begin the synergistic interaction critical to successful 
validation activities in their one-on-one discussions. Many felt that this 
workshop should be the forerunner of future workshops on specific topic 
areas because the workshop atmosphere provides a forum for discussing . both 
successes and failures which taken together often lead to more expeditious 
problem solutions. 
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PROGRAM 

Tues. . July 10. 1990 

8:30 Registration, Administration Building Auditorium 

9:00 Welcome Lonnie Reid 

9:10 Workshop Goals and Agenda Ray Gaugler 

9:20 Overview of NASA CFD Validation Program .... Ed Schairer 

9:45 LaRC Summary Scott Kjelgaard 

10:15 ARC Summary Joe Marvin 

10:45 Break 

11:00 MSFC Summary . . Paul McConnaughey 

11:45 LeRC Summary . Ray Gaugler 

12:15 Lunch - On your own, Lewis Cafeteria 

1:15 ARC Technical presentations .... William Lockman 

Paul McConnaughey 

2:15 LaRC Technical presentations Scott Kjelgaard 

3*15 Br63ik 

3:30 LeRC Technical presentations . . . Warren Hingst 

Michael Hathaway 

4:30 MSFC Technical presentation . . . Lisa Griffin 

5:00 Social Hour, Hors D'Oeuvres & Cash Bar, Ad. Building Foyer 

Wed.. July 11. 1990 

8:00 Instructions to working groups Ray Gaugler 

8:15 Meeting of working groups 

12:15 Lunch - On your own, Lewis Cafeteria 

1:15 Meeting of working groups 

3:30 Working groups report 

4:30 Adjourn 

Thurs.. July 12. 1990 

8:00 Government caucus 

12:00 Adjourn 
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2nd NASA CFD VALIDATION WORKSHOP NASA 

GOALS 

• SUMMARIZE PROGRESS SINCE '87 WORKSHOP 

• AFFIRM FUTURE DIRECTION OF NASA'S CFD VALIDATION PROGRAM 
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Langley Research Center 
CFD Code Validation 
Program Overview 


Scott O. Kjelgaard 

Experimental Methods Branch 
Fluid Mechanics Division 
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CFD Valid Workshop 


Outline 


» LaRC Approach to CFD Code Validation 

* Experimental/CFD Perceptions 

> CFD Code Validation Program Experiment 
Overview 

* Experiment Highlights 

* Concluding Remarks 
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CFD Code Validation 

Objective: Conduct detailed benchmark experiments yielding 
high-quality archival data bases. Use these data bases to 
validate CFD codes and develop empirical models of complex 
fluid dynamic phenomena 

Approach: 

Conduct cooperatively-designed experiments which provide: 




• detailed flowfield and surface flow measurements required 
by code developers 

• assessment of measurement errors 

• redundant measurements with more than one instrument, 
in more than one facility 

• documented results with archival storage of data for easy 
accessibility 


Acquire dedicated, advanced instrumentation systems for the 
primary code validation facilities 



Develop new instrumentation techniques when required 

July 1990 
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THE EXPERIMENTALIST'S VIEW? 
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CFD Code Validation Experiments 



LaRC Code Validation Program 

• 65° Swept Delta Wing 

• 76° Swept Delta Wing 

• Vertical Tail/Vortex Interaction Investigation 

• Rearward Facing Step Experiment 

• Particle Imaging Velocimeter (PIV) 

• Doppler Global Velocimeter (DGV) 

• Holocinematographic Velocimeter (HCV) 

• LTPT Laser Velocimeter 
Other Programs 

• F-18 Mean Flowfield Measurements 

• Transonic After-Body Experiment 

• Supersonic Coaxial Jet 

July 1990 
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76° Delta Wing Experiment 


Objective: 

Develop a detailed data base for the flow over sharped edged delta 
wings with a leading-edge sweep of 76°. Use three component LV 
to obtain flowfield measurements in the burst vortex over the 
delta wing. 


Approach: 

Conduct experiments in various facilities documenting the force 
and moments, surface pressures and flowfield velocities 
throughout the angle-of-attack range. 


Data acquired: 

•Forces, Moments 

•Static and Dynamic Surface Pressures 
Three component LV, 5-Hole Probe 




Facility 

BART, 7x10, Vigyan 

BART 

BART 
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Vertical Tail/Vortex Interaction 

Investigation 



Objective: 

• Provide a fundamental understanding of the 
vortex-fin interaction process 

• Determine the capability of CFD methods to predict 
the physics of the flowfield 

• Provide a data base for developing improved tail 
buffet criteria 

• Develop experimental measurement techniques 
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Vertical Tail/Vortex Interaction 

Investigation 


Annroach: 

• Cooperative MCAIR-Langley MOA in place to study 
vortex flows and vortex-fin interactions 

• Multi-disciplinary program involving structural 

dynamics and aerodynamics groups at MCAIR and 
Langley 

• Experiments to be conducted using a simple delta 
wing models with twin vertical tails 


July 1990 * 
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Vertical Tail/Vortex Interaction 





Experimental Set-up: 

• 76° delta wing using Hummel planform 

• 3 component LV measurements to determine the aerodynamic 
input forcing function to the vertical tails 

• Unsteady surface pressure measurements on the rigid tail to 
be correlated with the velocity fluctuations in the vortex 


Dynamic response from the flexible tail 


Measurement Techniques 

• Unsteady Surface Pressures 

• Tail Bending and Torsion 

• Laser Light Sheet Flow Visualization 

• 3 component Laser Velocimeter 
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65 ° 


Delta Wing Experiment 


Qbj^gtivgi 



Develop a detailed data base for the flow over rounded-edge 
delta wings with a leading-edge sweep of 65°. Measurements 
will be obtained throughout the Mach and Reynolds number 
range in LTPT, NTF, and 8' TPT at LaRC. 


Low Turbulence Pressure 
Tunnel 

• Static surface pressures 

• Flow visualization 

• 3 component LV 

8-foot Transonic Pressure Tunnel 

• Static surface pressures 

• Flow visualization 

• 3 component LV 
National Transonic Facility 

• Static surface pressures 

• Forces, moments, buffet 

• Vortex-core condensation 
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Geometric Features of NTF Delta Wing 

A te = 65°, c r = 2.14 ft, b = 2.0 ft 


Interchangeable leading Edges 

• [n e /c]s=.05,.15,.30% 

• 8 | x surface finish 


Sting 


• Symmetric 

• d/b = 0.14 


Flat plate center wing 

• t/c r = 0.034 

• 8|t surface finish 

• Sharp trailing edge 


Pressures: 


• 183 orifices / configuration 

• x/c r = .2, .4, .6,. 8,. 95 

• every 10% x! c r along l.e. 

• 0.010 inch diameter 
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Particle Imaging Velocimeter (PIV) 


)biective: 

Develop non-intrusive measurement techniques to obtain the 
data required for CFD code validation 


A pproach: 

Particle trajectories are tracked by double pulsing a laser 
light sheet and recording the images. Particles in the sheet 
form a double-spot pattern on the image when illuminated by 
the laser beam. The spacing and orientation of the spot-pairs 
provide an instantaneous measurement of the velocity in a 
two-dimensional plane in the flowfield 


July 1990 = 
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Photograph of the 
Particle Image Velocimetry 
Data Analysis 


July 1990 
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75° Delta Wing Investigation, a = 20.5° 

Flow Visualization 

• Surface using Ti02 

• Off-body using scanning laser light sheet 

Flowfield Surveys 

• Pitot pressure 

x/L = 0.3, 0.5, 0.7, 0.9, 1.1, Ee = 0.5, 1.0, 1.5 million 
Each survey plane contains ~ 2800 data locations 

• Velocity (5-Hole Probe) 

x/L = 0.7, 0.9, 1.1; Re = 0.5, 1.0, 1.5 million 
Each survey plane contains = 3400 data locations 

• Velocity (LV) 

x/L = 0.9; Re = 1.0 million 

Data obtained to evaluate measurement errors of 5-hole probe 


July 1990 
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COMPUTATIONAL GRID REFINEMENTS 

75° Delta Wing, x/L = 0.7 


65 X 65 (4225 Total) 



4 Patches (7222 Total) 



Single Dense Grid 


Embedded Grid 


wjA-l/^Ws 












STREAM WISE VELOCITY CONTOURS 

R e L = 500,000; a = 20.5°; x/L = 0.7 


Experiment 2 Level Refinement 



D88-1 & D88-4 Combined 


H3 






CFD Valid Workshop =— — — 

Concluding Remarks 



LaRC is dedicated to the development of a strong CFD Code 
Validation Program 

• Ongoing experimental program documenting various fluid 
dynamic phenomena 

• Experiments designed with experimental/computational 
cooperation 

• assessment of measurement errors 



• redundant measurements with more than one instrument, in more than 

one facility 

• documented results with archival storage of data for easy 

accessibility 


• Acquisition of the required instrumentation for CFD code 
validation facilities 

• Development of new instrumentation systems when required 


July 1990 
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OVERVIEW OF CFD FLOW MODELING AND VALIDATION 

AMES RESEARCH CENTER 


J. G. Marvin 


Second NASA CFD Validation Workshop 
Lewis Research Center 
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July 10-12, 1990 
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OUTLINE 


Background 

Experiments 

Highlights 

Instumentation 

Experiments 

Appendix 

Summaries And Key 
Results of Experiments 
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STAGES OF CODE DEVELOPMENT AND CORRESPONDING 
EXPERIMENTS RELATED TO PAYOFF IN AERONAUTICAL 
APPLICATIONS 


Science 

discipline 

Enabling 

technologies 

Research 

codes 

Pilot 

codes 

Production 

codes 

CFD 

• Algorithms 
•Grid generation 

• Computer Power 

•Technology 
integration 
• Limited 
pioneering 
applications 

• Wide range of 
applications 

• Applied in 
design 
environment 
•Cost effective 

EFD 

• Facilities 

• Instrumentation 

• Data acquisition 

Building block 
experiments 

Benchmark 

experiments 

Design 

experiments 

• Flow physics 
including 
numerical 
simulations 

• Flow modeling 

•Calibration 
• Validation 

• Configuration 
•Performance 
•System 
integration 


Payoff /J 
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KEY PROGRAM ELEMENTS FOR FLOW MODELING AND CFD VALIDATION 

AMES RESEARCH CENTER 


Focus On 3-D Reynolds Averaged Navier Stokes Code Validation 

Focus on Important CFD Applications 

Involve CFD'ers, Modelers And Experimentalists 

Design Experiments To Provide The Specific Data Necessary To Guide 
And Verify CFD Computations 


FACILITY AND INSTRUMENTATION ENHANCEMENTS 

1986-PRESENT 


FACILITIES REACTIVATED FOR HYPERSONICS 
3.5' Hypersonic Wind Tunnel 
Ballistic Range 

Combustion Driven Shock Tunnel 
Electric Arc Shock Tube 

INSTRUMENTATION ADVANCES/DEVELOPMENTS 
New Near-Wall measurement capability using LDV 
New LDV Systems for Hi Re Channel And 3.5' HWT 
Pressure Sensitive Chemical Luminescence Paints 
LIF For Temperature And density Measurements in Wind Tunnelsu 
Holographic System For Ballistic Range Flow Field Density Measurements 
LIF For Combustion Driven Shock Tunnel-Gas Species, Tempertures 


FLOW MODELING AND CFD VALIDATION EXPERIMENTS 

1990 


BUILDING BLOCK EXPERIMENTS BENCHMARK EXPERIMENTS 


EXPERIMENT 

CONTACT 

ISSUE 

EXPERIMENT 

CONTACT 

CODE 

1) Adv. Press.Grad 

Watmuff 

Flow Physics 

1 2) Low Aspect Ratio 

Olsen 

TNS 

Boundary Layer 



Wing 



2) Rearward Step 

^Jovic 

Flow Physics 

1 3) Airfoil w/wo Pas- 

Mateer 

TURF 



sive Shock Control 



3) 3-d Spin Flow 

Driver 

Turb. Model 



- 



14) Generic All Body* 

Lockman 

UPNS 

4) 3-d Wedge Flow 

Johnson 

Turb. Model 





15) Ballistic Range- 

Strawa 

TUFF/STUFF 

5) Wing Tip Vortex 

Zilliac 

Turb. Model 

Cone Real Gas 



6) 3-d supersonic 

Horstman 

Turb. Model 

16) NASP Nozzle 

Lockman 

F3D 

Shock Interaction 



17) STOVL 

Van Dalsem 

F3D 

7) SupersonicShock 

Settles 

Flow Model 




Vortex Interaction 

(Penn St.) 


18) Transonic WT 

Roberts 

TNS 



Wall Boundary 

(Stanford) 


8) SSME Turn-Around 

Monson 

Turb. Model 

Conditions 



Duct* 






9) Coanda 

Brown 

Turb. Model 




10) Hypersonic Shock 

Horstman 

TUrb. Model 




Interaction 






11) Shock Tube 

Sharma 

Reaction Rates 





* Chosen For Technical Presentation 
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KEY ACCOMPLISHMENTS 
Instrumentation 
Experiments 
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LASER DOPPLER VELOCIMETER FOR 3-D NEAR-WALL MEASUREMENTS 



To photodetector 


Run 

O 511,4/90 
A 491,4/90 

□ 471,4/90 — Spalart 

V 121,3/90 ------ Kartsson and Johansson 

O 181, 4/90 ' — 1 Barlow and Johnston 
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LDV FEASIBILITY STUDY IN A HYPERSONIC SHOCK-WAVE/I^^ 
' TURBULENT-BOUNDARY-LAYER INTERACTION ' 


■ : ‘ -t^-p-v '. ■ $ «5^« | | ~V' /■ 

. - -i-' - . v- ^ ^ V --.-- '■' 
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LASER-INDUCED FLUORESCENCE (LIF) 

MEASURES TEMPERATURE, DENSITY, SPECIES 




EXPERIMENTAL 

CONFIGURATION 


MOLECULAR ENERGY LEVELS 
AND RADIATIVE FREQUENCIES 


5S 
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ENERGY/FREQUENCY 


TEMPE 

CONST/I 


S 1 * 

(0 

E 

CB 

5 

« 1.2 
c 

<D 

TJ 


0 









KEY ACCOMPLISHMENTS 


Instrumentation 



Experiments 
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I VI IVDt\ NAME 




ORIGINAL PAGE SS 
OF POOR QUALITY 




TOPOLOGY OF COHERENT STRUCTURES 

Reg = 2000 Uqo = 10 m/sec 


Y-RAKE 



U c t/5 
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ENSEMBLE AVERAGE 

X/h = -2.6 










VERIFICATION EXPERIMENT FOR TRANSONIC FLOW CODES 


SUPERCRITICAL AIRFOIL WITH LOV 






'P' ; 




Vi. 


■m 


C 


ttTTV 


w 










■ i 

A 


€ 


Vk7. f V- 

!|>i 

3* 


• VrM 


m 


TEST CONDITIONS 

‘■.V# 0 . 73 CM < 4*0 

2 X 10 ® < Re < 6 X 10 6 
0 . 5 <a< 1 £ 

MEASUREMENTS 

’ PRESSURE Ip. <P1) _ 

FLOW FIEU) IU.;V. u‘ V . . .) 

^ OIL FLOWt 

V LIFT. DRAG d 


WING PRESSURES 
M - 078 « - 1° 



EXPERIMENT 

ADVANCED TURBULENCE MODEL 
STANDARD TURBULENCE MODEL 


WAKE VELOCITY 
■ M - 078 a -*'7° Kff 1.04 
* ■ , J' „ l s;* 

O EXPERIMENT' ~ 

— ADVANCED TURBULENCE 
MODEL 

ST ANDARO TURBULENCE 

MODEL 
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RANS TRANSONIC AIRFOIL PREDICTIONS 

TURBULENCE MODELING IMPROVEMENTS RAE 2822 AIRFOIL 


□ EXPERIMENT 

NONEQULIBRIUM MODEL 

EQUILIBRIUM MODEL 


ATTACHED CASE 
M « 0.73, a«3° 


SEPARATED CASE 
I « 0.75, a«3° 


-Aft 8 


.4 .6 

X/C 


1.0 0 


.4 6 

X/C 




nn M*rm m**m (raw 







EFFECT OF TURBULENCE MODEL ON WII 
ONER A M-6 WING, M*0.84 o 


PRESSURE CONTOURS 



BALDWIN-LOMAX 


oo EXPERIMENT 



X/C 


- 1.6 

- 1.2 


-.8 


-.4 

0 


.4 


.8 

1.2 




PRESSURES 




oo EXPERIMENT 



0 .2 .4 .6 .8 1.0 


x/c 
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Pitching 

moment coefficient 


REAL-GAS FORCE AND MOMENT MEASUREMENTS AND CFD 
CORRELATION FROM THE BALLISTIC RANGE 

Re = 10®, V^s 5 km/sec, Mn = 14.5 


.2 


Measurement 

c 

<D 

O 1 

1 

(Error: ± 3%) 

1 .1 



u 

O 



O) 

£ 


Computations { 

6 




1 


Non-equilibrium air 

I u 


2.5 5.0 7.5 

a , Angle-of-attack, deg 



1 2 3 4 5 6 7 8 

a , Angle-of-attack, deg 


g 10 

o 

Q. 

S c 

CO 3 

■o 


o 

z 




• ••••• 

1 I T T t ! i 


1 2 3 4 5 6 7 

(X , Angle-of-attack, deg 


Shadowgraph of blunt cone 



Computed pressure contours 



* 


b 7 
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National Aeronautics and 
Space Administration 


r 


X--609 
His 23 


NASA - Ames Research Center 
Hypervelocity Free Flight Facility - March 1990 


FINITE FRINGE INTER FRROGR AM 


5 degree half-angle blunt cone 



NASP Nozzle/Afterbody Computations 

Experimental Geometry 3-D Composite Grid Topology ; i 1|| 


1h ; :»i . ^ 
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Mach Contours in Symmetry & Outflow Planes 


Pressure Contours & Particle Traces 

, S.M. Ruffin 8c E. VenRatapathy 
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STATUS OF MSFC CFD APPLICATION AND VALIDATION ACTIVITIES 


BY 


L.A. SCHUTZENHOFER, H.V. McCONNAUGHEY, P.K. McCONNAUGHEY 


COMPUTATIONAL FLUID DYNAMICS BRANCH 
AEROPHYSICS DIVISION 
STRUCTURES AND DYNAMICS LABORATORY 
NASA/MARSHALL SPACE FLIGHT CENTER 


i ■ 

i 


SECOND NASA CFD VALIDATION WORKSHOP 
JULY 10-12, 1990 
NASA LEWIS RESEARCH CENTER 


STATUS OF MSFC CFD APPLICATION AND VALIDATION ACTIVITIES 

OVERVIEW 

• FOCUS OF MSFC CFD ACTIVITIES 

• DESIGN APPLICATIONS t 

• VALIDATION ACTIVITIES 

• CFD REQUIREMENTS 

... 

• CFD EXPECTATIONS 


•SUMMARY 


STATUS OF MSFC CFD APPLICATION AND VALIDATION ACTIVITIES 

FOCUS OF MSFC CFD ACTIVITIES 


• SUPPORT PROGRAM OFFICES 

- " QUICK TURNAROUND " APPLICATIONS; SENSITIVITY ANALYSIS ; RETROFITTABLE D E SI GN OPTIONS 

- INTERACT WITH HARDWARE CONTRACTORS IN DEVELOPMENT OF DESIGN ENVIRONMENTS 

- PROVIDE " SMART BUYER " CAPABILITY FOR LONG-TERM APPLICATIONS 

- DEVELOP SUBSYSTEMS CFD MODELS 

- FOCUS MSFC CFD ACTIVITIES/ PROVIDE CENTERWIDE CFD SUPPORT 

• FOCUS DEVELOPMENT OF CFD METHODOLOGY 

- INTERACT WITH ARC. LeRC. LaRC. AND OTHER RESEARCH ORGANIZATIONS IQJFQC11S 
TECHNOLOGY DEVELOPMENT TOWARDS MSFC HARDWARE RELATED PROBLEMS 

- DEVELOP REQUIREMENTS FOR CFD CODE VERIFICATION 

- VERIFY CODES THROUGH BENCHMARK COMPARISONS 
-ADVANCE CFD TECHNOLOGY FOR APPLICATIONS 

• DEVELOP ADVANCED HARDWARE TECHNOLOGY CONCEPTS 

-TURBINE STAGE 

- PUMP STAGE 

- NOZZLES, PREBURNERS, ETC. 
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STATUS OF MSFC CFD APPLICATION AND VALIDATION ACTIVITIES 

DESIGN APPLICATIONS 

• PROGRAM SUPPORTING ACCOMPLISHMENTS 

• DESIGN APPLICATION EXAMPLES 

— atd discharge volute side load analysis 

—MAIN INJECTOR LOX INLET FLOW INDUCED VIBRATION 
—NOZZLE/MAIN COMBUSTION CHAMBER (MCC) MISMATCH 


• CURRENT PROGRAMS 


STATUS OF MSFC.CFD APPLICATION AND VALIDATION ACTIVITIES 

PROGRAM SUPPORTING ACCOMPLISHMENTS 

INHQUSE CONT. 


SSME • HPFTP TURBINE BLADES X 

• TURBINE DISK CAVITIES X 

• LOX PUMP BEARING INLET CAVITY X 

• LOX PUMP BEARINGS 

• FUEL PREBURNER X 

• LOX PREBURNER X 

• LOX MANIFOLD TEE (4000 Hz) X 

• HOT GAS MANIFOLD/MANIFOLD STRUTS X 

• PUMP COOLANT FLOW PATHS X 

• NOZZLE/MCC MISMATCH X 

• HPOTP NOZZLE PLUG TRAJECTORIES X 

• TRANSIENT BEHAVIOR OF FUEL PREBURNER MANIFOLD 

• BEARING DEFLECTOMETER (TTBE) X 

• ENGINE 0212 BEARING CAGE DEBRIS X 

• ENGINE 0209 MCC COOLANT LINER CRACKING X 

• ENGINE 2019 HPFTP STATIC SEAL PARTICLE 

TRAJECTORIES X 

ATD • TURBINE INLET TEMP. REDISTRIBUTION X 


TURBINE TEMP. PROFILE REDISTRIBUTION 

ROTOR-STATOR INTERACTION 

TURNAROUND DUCT AND HOT GAS MANIFOLD 

BEARING ANALYSIS 

LOX PUMP INLET SCROLL 

FUEL PUMP INTERSTAGE CROSSOVER DUCTS 

FUEL PUMP INLET SCROLL 


• LOX PUMP DISCHARGE VOLUTE X 

• SEALS 

• HOT GAS MANIFOLD CHECKOUT CHAMBER X 

• FUEL PUMP DISCHARGE VOLUTE X 

• HPFTP AFT DISK CAVITY X 

• HPOTP TURNAROUND DUCT VANE FAILURE X 

• UNSTEADY MULTISTAGE TURBINE LOADS X 




Mm 









X 


X 

X 

X 

X 

X 


X 

X 

X 

X 

X 

X 

X 

X 

X 

X 
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STATUS OF MSFC CFD APPLICATION AND VALIDATION ACTIVITIES 

PROGRAM SUPPORTING ACCOMPLISHMENTS 


SRB 


AFE 

SPACE STATION 

ADVANCED PROGRAM 
DEVELOPMENT 


INHO USE 


C ONT . 


BORE FLOW 

— CANTED NOZZLE 

— BROKEN INHIBITOR 

FIELD JOINT 

— FLOW AND THERMAL TRANSIENT 

— PRESSURIZATION TRANSIENT 

NOZZLE-TO-CASE JOINT 

— FLOW AND THERMAL TRANSIENT 

— PRESSURIZATION TRANSIENT 

MNASA MOTOR 

— WITH BLAST TUBE 

— WITHOUT BLAST TUBE 


X 

X 

X 

X 

X 

X 


X 

X 

X 

X 

X 

X 


• AEROTHERMAL ENVIRONMENTS 
— DSMC 

— NS 

• CONTAMINATION TRACKING 

• ECLSS MODULE AND NODE 


X 

X 

X 

X 


• EXTERNAL TANK GAMMA RAY IMAGING 
TELESCOPE 
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VISCOUS ANALYSIS OF THE ATD OXYGEN AND FUEL 
PUMP DISCHARGE VOLUTES 


• OBJECTIVE: TO OBTAIN SIDELOADS DUE TO HIGH PRESSURE PUMP VOLUTES. 

• JUSTIFICATION: ACCURATE SIDELOAD PREDICTION IS CRITICAL FOR SUCCESSFUL BEARING 
DESIGN. WATER-FLOW RIG DATA AVAILABLE FOR OXYGEN PUMP VOLUTE. NO DATA 
AVAILABLE FOR FUEL PUMP VOLUTE UNTIL 3RD HOTFIRE TEST SERIES OF TURBOPUMP. 

• APPROACH: USED FDNS2D TO PERFORM INCOMPRESSIBLE, 2-D TURBULENT ANALYSIS OF 
VOLUTE GEOMETRY. JACOBIAN MODIFIED TO ACCOUNT FOR 3-D CROSS-SECTION AREA. 

USING 511 X 97 GRID POINTS FOR 5-VANED, DUAL EXIT LOX VOLUTE. USING 481 X 79 GRID FOR 
1 3-VANED SINGLE EXIT FUEL VOLUTE. 

• RESULTS/IMPACT: 2-D RESULTS: LOX VOLUTE ANALYSES MATCHED RESULTANT LOAD (700 
LBS.) BUT NOT DIRECTION. FUEL VOLUTE RESULTS INDICATE VERY SMALL LOADS (< 100 LBS.). 
PSEUDO 3-D ANALYSIS STILL UNDERWAY. 
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STATUS OF MSFC CFD APPLICATION AND VALIDATION ACTIVITIES 

MAIN INJECTOR LOX INLET CFD RESULTS 


FLIGHT CONFIGURATION Re = 10 6 ; St = .375 (3800 Hz) 



TIME (TxUREF/XREF) 

B PRESSURE/LOAD VARIATION WITH TIME 







STATUS OF MSFC CFD APPLICATION AND VALIDATION ACTIVITIES 


NOZZLE/MCC MISMATCH FLOW FIELD ANALYSIS 


OBJECTIVES 

• DETERMINE TEMPERATURE IN G1 5 CAVITY FOR VARIOUS MISMATCH CONFIGURATIONS 

• ASSESS SENSITIVITY OF CAVITY TEMPERATURE TO VARIATIONS IN MAXIMUM 
PROTRUSION, CIRCUMFERENTIAL DIFFERENTIALS IN PROTRUSION, AND GAP 

• DEVELOP ACCEPT/REJECT CRITERIA FOR FUTURE ENGINE CONFIGURATIONS 


APPROACH 

• CFD SIMULATIONS OF 2D, CLOSED CAVITY FLOW (APPROXIMATE GEOMETRY) 

- ASSESS CAVITY PRESSURE AS A FUNCTION OF PROTRUSION 

- DEDUCE APPROXIMATE CIRCUMFERENTIAL FLOW AS A FUNCTION 
OF MAXIMUM PROTRUSION 

• CFD SIMULATIONS OF 2D CAVITY WITH MASS FLOW OUT OF CAVITY BASE 
(EXACT GEOMETRY) 

- PERFORM PARAMETER SENSITIVITY STUDY 

- DETERMINE FOR WHAT PARAMETER VALUES THE CAVITY TEMPERATURE 
REACHES THE CRITICAL TEMPERATURE AT WHICH HEAT CONDUCTION 
AT NOZZLE LIP WILL CAUSE BLUEING OF SEAL 


PROGRESS IMPACT 

• RESULTS COMMUNICATED TO CHIEF ENGINEER FOR SSME AND TO HARDWARE 
CONTRACTOR 

• RESULTS USED IN DEVELOPMENT OF ACCEPT/REJECT CRITERIA 


1 !! 
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STATUS OF MSFC CFD APPLICATION AND VALIDATION ACTIVITIES 

CURRENT PROGRAMS 


] 


Development Programs 


1 KEY PROGRAM MILESTONES 
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Advanced Launch System 
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Advanced Solid Rocket Motor 
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CFD ACTIVITIES 



















Aeroassist Fliaht Experiment 
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Space Station Freedom 
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ECLSS 
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Alternate Turbopump 
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VALIDATION ACTIVITIES 


• BENCHMARK VALIDATION PLAN 

• LEVEL I: FUNDAMENTAL FLOWS 

• LEVEL II: SUBCOMPONENTS 

• LEVEL III: INTERACTIVE COMPONENTS/SYSTEMS 

• END PRODUCTS 


• CONSORTIUM 




STATUS OF MSFC CFD APPLICATION 
VALIDATION ACTIVITIES 


• LEVEL I: FUNDAMENTAL FLOWS 


FLAT 

CHANNEL 

STRAIGHT 

CODE PLATE 

FLOW 

PIPE 

INS3D 

1 

1 

INS3DLU 1 

2 

1 

INS3DUP 1 

1 

2 

CMINT 

2 

1 

FDNS 3 

3 

4 

REFLEQS 1 

3 

1 

MAST 


• - 

* NUMBER OF PEOPLE WORKING PROBLEM 




CIRCULAR 

RECT 

BACKWARD 

PIPE BEND 

PIPE BEND 

✓ 

FACING STEP 

1 
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2 

1 

2 

2 

2 

2 

2 

1 

2 

4 

2 
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2 

1 
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STATUS OF MSFC CFD APPLICATION AND VALIDATION ACTIVITIES 

BENCHMARK VALIDATION PLAN 


• LEVEL 2 - SUBCOMPONENTS 

CASCADE 

DISK CAVITY 

SINGLE INJECTOR 

UNSTEADY FLOW OVER CYLINDER 

AXISYMMETRIC DIFFUSER 

CURVED DUCT WITH TURNING VANES 

HIGH CURVATURE DUCT WITH SE PE RATION 

ROTATING DUCT 

ROTATING CURVED DUCT 

NOZZLE 

CYLINDER WITH BLOWING WALLS 
SLOTTED CYLINDER WITH BLOWING WALLS 


• LEVEL 3 - INTERACTIVE COMPONENTS 
OR SYSTEMS 

HOT GAS MANIFOLD PILOT MODEL 
UNSTEADY TURBINE STAGE 
MAIN INJECTOR/NOZZLE/PLUME 
UNSTEADY PUMP STAGE 
INLET AND DISCHARGE VOLUTES 
AEROASSIST FLIGHT EXPERIMENT 
BEARING CAVITY MODEL 


STATUS OF MSFC CFD APPLICATION AND VALIDATION ACTIVITIES 

BENCHMARK VALIDATION PLAN 


• DOCUMENTATION REQUIREMENTS 


• OBJECTIVE: 


DESCRIBE PROBLEM AND THE FOCUS/EXPECTATION OF 
THE CFD COMPUTATION 


• COMPUTATIONAL PROBLEM DEFINITION, PROBLEM GOAL, BOUNDARY 
STRATEGY: CONDITIONS, INITIALIZATION, USER ISSUES 


• DATA BASE: 


DEFINITION OF GEOMETRY, Re, M, DATA SOURCE, ETC. 


• COMPUTER CPU TIME, STORAGE REQUIREMENTS, SSD REQUIREMENTS 

REQUIREMENTS: 


• LESSONS LEARNED: DESCRIPTION OF FALSE STARTS, GRID SENSITIVITY, 

INITIALIZATION, B.C. PROBLEMS. 

• SUMMARY: STRENGTH/WEAKNESS IN COMPUTATION VS. DATA 

COMPARISON. EVALUATION OF CODE FOR POTENTIAL 
HARDWARE APPLICATION. 


DOCUMENTATION SHOULD BE CONCISE . 


STATUS OF MSFC CFD APPLICATION AND VALIDATION ACTIVITIES 

BENCHMARK VALIDATION PLAN 

• END PRODUCTS: 

- DIMENSIONLESS EVALUATION CRITERIA 

- ERROR ESTIMATION METHODOLOGY 

- ASSESSMENT OF CODES BASED UPON CRITERIA AND ERROR ESTIMATION 

• FOCUS OF CRITERIA: 

- FLOW RANGES (MACH NUMBER, REYNOLDS NUMBER ~) 

-CODE (SOLUTION ALGORITHM, CODE ARCHITECTURE, BOUNDARY CONDITIONS ~) 

- FLOW PROCESS MODELS (TURBULENCE MODEL, MULTISPECIES --) 

- USER SPECIFIC EFFECTS (GRID GENERATING, INITIALIZATION ~) 
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STATUS OF MSFC CFD APPLICATION AND VALIDATION ACTIVITIES 

CONSORTIUM 


• CONSORTIUM OBJECTIVES 

• FOCUS CFD APPLICATIONS IN PROPULSION 

•TECHNOLOGY ACQUISITION PHASE 

•DIRECT BASELINE PROGRAM TOWARDS IMPROVED ACCURACY, 

STABILITY, AND EFFICIENCY 

•LARGE SCALE SUBSYSTEM TECHNOLOGY VALIDATION 

•STIMULATE CFD VALIDATION TOWARDS PROPULSION FLOWS 

•DIRECT APPLICATIONS CODES TOWARD DESIGN TOOLS AND 
ADVANCED TECHNOLOGY HARDWARE CONCEPTS 

• IDENTIFY NATIONAL CFD PROPULSION REQUIREMENTS 

• STIMULATE A FORUM FOR GOVERNMENT, INDUSTRY, AND UNIVERSITY 
INTERACTIONS 

• ENCOURAGE INDUSTRY TO PARTICIPATE IN CFD DEVELOPMENT WITH IRAD FUNDS 

• PROVIDE SYNERGISM IN THE CFD COMMUNITY 

• PROVIDE PEER REVIEW OF CFD PROGRAMS 


STATUS OF MSFC CFD APPLICATION AND VALIDATION ACTIVITIES 


i CONSORTIUM FOR CFD APPLICATIONS IN PROPULSION TECHNOLOGY 

ORGANIZATIONAL STRUCTURE 
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STATUS OF MSFC CFD APPLICATION AND VALIDATION ACTIVITIES 

ADVANCED TECHNOLOGY HARDWARE 

• CONSORTIUM TEAMS 

FOCUS DEVELOPMENT OF CFD METHODOLOGY 

AND 

DEVELOP ADVANCED TECHNOLOGY HARDWARE CONCEPTS 
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STATUS OF MSFC CFD APPLICATION AND VALIDATION ACTIVITIES 


• CONSORTIUM TEAM VALIDATION ACTIVITIES 


TURBINE TEAM 

- CASCADES 

- ROTOR STATOR 

- SSME 

PUMP TEAM 

- ROCKETDYNE INDUCER 

COMPLEX FLOW PATHS 

- DISK CAVITIES 

• HOM PILOT MODEL 

- DUCT FLOWS 


- SRM TEAM 

- CSD CYLINDER 

- SLOTTED CYLINDER 

- SRM PILOT MODEL 

• COMBUSTION TEAM 
■ H .MIXING 

• DUMP COMBUSTION 

• NOZZLE FLOW 


STATUS OF MSFC CFD APPLICATION AND VALIDATION ACTIVITIES 

CONSORTIUM 


• TURBINE STAGE TECHNOLOGY TEAM 


OBJECTIVE- DEVELOPMENT, ENHANCEMENT, VALIDATION, AND DEMONSTRATION OF CFD 

TOOLS FOR TURBINE STAGE DESIGN 

APPROACH: • ENHANCE EXISTING DATA BASE (COLD FLOW RIG AND HOT-FIRE SYSTEM 

TESTS), DEVELOP/ENHANCE CFD MODELS AND CODES, 

EXPERIMENTALLY DEMONSTRATE CFD-PREDICTED DESIGN 
IMPROVEMENTS 

• ASSEMBLE TURBINE TECHNOLOGY EXPERTS REPRESENTING DIFFERENT 
AREAS OF EXPERTISE. POOL TECHNICAL RESOURCES AND EXPERIENCE 
BASES TO BRING VARIOUS PESPECTIVES AND ABILITIES TO BEAR ON 
TURBINE STAGE ISSUES AND ON ACCOMPLISHMENT OF TEAM 
OBJECTIVE 

• DEVELOP AND IMPLEMENT PLAN WHICH COORDINATES 
MSFC-SUPPORTED TURBINE STAGE ANALYSIS AND TEST ACTIVITIES 

• IDENTIFY MILESTONES, DELIVERABLE PRODUCTS, AND TARGET DATES 

• IDENTIFY ADDITIONAL ACTIVITIES REQUIRED TO ACCOMPLISH TEAM 
OBJECTIVE 

• INTERACT THROUGH REGULAR MEETINGS TO STATUS, CROSS-FERTILIZE, 
CRITIQUE, AND/OR DIRECT EACH ACTIVITY 


STATUS OF MSFC CFD APPLICATION AND VALIDATION ACTIVITIES 

CONSORTIUM 


• TURBINE STAGE TECHNOLOGY TEAM 

STATUS: 5 MEETINGS HELD TO DATE-4/28/89, 8/23/89, 1 1/1/89, 1/17/90, 

4/19/90 


MEMBERSHIP: 


• NASA MARSHALL SPACE FLIGHT CENTER (MSFC) 

• NASA AMES RESEARCH CENTER (ARC) 

• NASA LEWIS RESEARCH CENTER (LeRC) 

• ROCKETDYNE (RDYN) 

• PRATT 8, WHITNEY (P&W) 

• AEROJET 

• GENERAL ELECTRIC COMPANY 

• UNITED TECHNOLOGIES RESEARCH CENTER (UTRC) 

• CALSPAN— UNIVERSITY OF BUFFALO RESEARCH CENTER 
(CUBRC) 

• ROTODATA, INC. 

• SCIENTIFIC RESEARCH ASSOCIATES (SRA) 

• UNIVERSITY OF ALABAMA— HUNTSVILLE (UAH) 

• PENSYLVANIA STATE UNIVERSITY (PSU) 

• MISSISSIPPI STATE UNIVERSITY (MSU) 


STATUS OF MSFC CFD APPLICATION AND VALIDATION ACTIVITIES 


CONSORTIUM - TURBINE TEAM 
KEY MILESTONES 



KEY* 1 

2. Validation against unsteady aerodynamic data 

3. Validation against SSME or advanced concept data 

4. Steady and unsteady aerodynamic data, unsteady heat transler data 

5. ATD steady aerodynamic data 


5-11344-0-310 









STATUS OF MSFC CFD APPLICATION AND VALIDATION ACTIVITIES 

CONSORTIUM - TURBINE TEAM 

• BASELINE TURBINE MIDSPAN CONTOURS/CALCULATED STREAMLINES 


1st VANE 1st BLADE 2nd VANE 2nd BLADE 



NOTE: ROW TO ROW SCALES ARE DIFFERENT 






STATUS OF MSFC CFD APPLICATION AND VALIDATION ACTIVITIES 

CFD REQUIREMENTS 

• GRID GENERATION 

- ELECTRONIC SURFACE GENERATORS COMPATIBLE WITH CFD CODES 

- INTERACTIVE GRID PACKAGES WITH INTERROGATION SCHEMES 

- ADAPTIVE GRID PACKAGES 

• EVALUATION CRITERIA AND ERROR ANALYSIS METHODOLOGY 

- STANDARDIZED DIMENSIONLESS PARAMETERS 

• VALIDATION EXPERIMENTS • £' 

•CODES 7 

- BENCHMARKING USING STANDARDIZED METHODS 

- MULTIBLOCKING OR ZONAL STRATEGY 

• FLOW PROCESS MODELS 

-TURBULENCE AND TRANSITION 

- MULTISPECIES 

- CHEMISTRY WITH AND WITHOUT REACTIONS 

- MULTIPHASE 
-ATOMIZATION 

- EVAPORATION 
-MIXING 

• POSTPROCESSING 

- DEMONSTRATION OF CONSERVATION PROPERTIES 

- ERROR ANALYSIS 

• GRAPHICS 

-INTERACTIVE 

-THREE DIMENSIONAL FLOW DESCRIPTIONS 
-ANIMATION 

5 - 9807 - 0-238 
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STATUS OF MSFC CFD APPLICATION AND VALIDATION ACTIVITIES 

CFD EXPECTATIONS 


• DIRECT HARDWARE DESIGN UTILIZING CFD 

— PROVIDE INITIAL IMPACT IN DESIGN 

— PERFORM DESIGN OPTIMIZATION STUDIES 

— DEVELOP ADVANCED TECHNOLOGY HARDWARE CONCEPTS 


• BENCHMARKED/VALIDATED CODES 

— LAMINAR FLOWS 

— TURBULENT FLOWS ut .p 

— ACOUSTIC PROBLEMS 

— CERTAIN CLASS OF UNSTEADY PROBLEMS 

• GENERAL METHODOLOGY FOR ERROR ESTIMATION AND EVALUATION CRITERIA 
FOR CODES 

• USER FRIENDLY CODES 

— B.S. LEVEL ENGINEER 2-3 YRS EXPERIENCE 

— GUIDELINES FOR CLASSES OF PROBLEMS 

— CAD/CAM/CAE; GEOMETRY GRID GENERATION 

— GENERALIZED BOUNDARY CONDITIONS 

— ALGORITHM/GRID OPTIMIZATION FOR SOLUTION EFFICIENCY 

— ARTIFICIAL INTELLIGENCE/EXPERT SYSTEMS 

— MODULAR CODES 

• FLOW ADAPTIVE GRIDS FOR CURRENT CLASS OF PROBLEMS 

• MULTIPLE SCALE AND/OR ZONAL TURBULENCE MODELS, MULTIPHASE, 
MULTISPECIES, COMBUSTION FLOW PROCESS ENGINEERING MODELS 
EVOLVED FROM EXPERIMENTS AND CFD ANALYSIS 


STATUS OF MSFC CFD APPLICATION AND VALIDATION ACTIVITIES 

SUMMARY 


• CFD BEING APPLIED IN THE HARDWARE DESIGN PROCESS AT MSFC 

• FOCUS ON TEAM STRATEGY; RIGHT PROBLEM 

• QUICK TURN-AROUND; SENSITIVITY ANALYSIS 

• RETROFITTABLE DESIGN OPTIONS 


• VALIDATION ACTIVITIES PROVIDING COMMUNITY SYNERGISM 

• LEVEL I BENCHMARK PLAN IN PLACE / IN PROGRESS 

• CONSORTIUM TEAMS INTEGRATE ADVANCED HARDWARE DESIGN AND 
VALIDATION ACTIVITIES 

• FOCUS CFD GOALS AND DEFINE REQUIREMENTS 

• INTEGRATE ACTIVITIES, POOL RESOURCES, AND TRANSFER 
TECHNOLOGY BETWEEN GOVERNMENT, CONTRACTORS/AND UNIVERSITIES 

• CFD EXPECTATIONS ARE APPLICATIONS-ORIENTED AND REQUIRE SUSTAINED 
TECHNOLOGY TRANSFER FROM RESEARCH CENTERS 
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National Aeronautics and 
Space Administration 


2nd NASA CFD 


VALIDATION WORKSHOP 



NASA 


Lewis Research Center 



SUMMARY OF LEWIS RESEARCH CENTER VALIDATION EFFORTS 


Raymond E. Gaugler 

Chief, Turbomachinery Flow Physics Branch 


National Aeronautics and 
Space Administration 

Lewis Research Center 


CLOSELY COUPLED EXPERIMENTAL 
AND COMPUTATIONAL RESEARCH 



I 


UNDERSTANDING OF FLOW PHYSICS 
ACCURATE PREDICTIVE COOES 


to 1 



National Aeronautics and 
Space Administration 

Lewis Research Center 


INTERNAL FLUID MECHANICS DIVISION 
HEAT TRANSFER BRANCH 


NASA 


LINEAR TRANSONIC CASCADE 








National Aeronautics and 
Space Administration 


INTERNAL FLUID MECHANICS DIVISION 

i ii— at — rn AMOCCD DDAMPU 


NASA 



LINEAR TRANSONIC CASCADE 


INLET BOARD 
POSITION INDICATOR 



INLET 

BOARDS 


INLET 

BOARD 

DRIVE 


DISC 

DRIVE 


J 



NASA 

C-87-09769 
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HEAT TRANSFER COEFFICIENTS 
ALONG DUCT CENTERLINE 



EXPERIMENTAL AND PREDICTED HEAT TRANSFER COEFFICIENTS 


CD-87-30'3 U 


NASA 

C-84-6765 





LOW SPEED CENTRIFUGAL COMPRESSOR 

BENCHMARK EXPERIMENTS FOR COMPUTATIONAL CODES 




OF POOR QUALITY 




Notional Aeronautics and 
Spoce Administration 

Lewis Research Center 


CFD VALIDATION PROGRAM 




TURBOMACHINERY BLADE ROW INTERACTIONS 

£ ■ ■ 

MULTI-STAGE COMPRESSOR FLOW PHYSICS 


ill 


Notional Aeronautics and 
Spoce Administration 

Lewis Research Center 




MULTISTAGE COMPRESSOR FLOW PHYSICS PROGRAM 













National Aeronautics and 
Space Administration 

Lewis Research Center 


CFD VALIDATION PROGRAM 





TURBOMACHINERY BLADE ROW INTERACTIONS 


MOTIVATION - 

Current Designs Based On Time-Averaged 
Axisymmetric Analyses 


Actual Flow Field Highly Non— Axisymmetric 
And Time Dependent 


OBJECTIVE 


Experimentally Validated Numerical Analysis 
Of Flows In Multi-Stage Turbomachinery 

Improved Understanding of Flow Physics 
Of Blade— Row Interactions 
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Notionol Aeronautics and 
Space Administration 

Lewis Research Center 


NASA 


LOW SPEED MULTISTAGE AXIAL-FLOW COMPRESSOR 

BENCHMARK EXPERIMENTS FOR COMPUTATIONAL CODES 


ft. diam. AXIAL-FLOW COMPRESSOR ASSEMBLY 





National Aeronautics and 
Space Administration 

Lewis Research Center 



Notional Aeronautics and 
Space Administration 

Lewis Research Center 


NASA 



M ULTISTAGE C O MPRESSOR FLOW PHYSICS PROGRAM 

SUMMARY 


A SYNERGISTIC PROGRAM PULLING .TOGETHER CFD EXPERIMENTS 

capabiuty SPEED interactive computing to provide a "Sd 


• CFD - UTILIZE AVERAGE-PASSAGE METHOD FOR MULTI-STAGE 

SIMULATION AND TIME ACCURATE METHODS TO SUPPORT 
CLOSURE MODELING 


E XPERIMENTS - INITIAL RESULTS FROM HIGH-SPEED, SINGLE 

STAGES - CL OSURE MODELING AND VALIDATION/ 
CALIBRATION FROM EXISTING HIGH-SPEED 
MULTISTAGE AND NEWLY REQUIRED LOW SPEED 
MULTISTAGE 

• M EASUREMENT TECH. - INSTRUMENTATION ADVANCES REQUIRED 

FOR 3D VELOCITY AND SCALAR 
MEASUREMENTS IN COMPRESSIBLE FLOWS 

H IGH-SPEED INTERACTIVE COMPUTING - REQUIRED FOR TIMELY CFD 

AND EXPERIMENTAL 
INTERACTION 


LeRC/lFMD INLET DUCT AND NOZZLE HIGH SPEED 

VALIDATION EXPERIMENTS 

• VALIDATION WORKSHOP EXPERIMENTS 

- CROSSING SHOCKS/BOUNDARY-LAYER INTERACTION 

- LeRC/UNIVERSITY UNSTEADY SHOCK/BOUNDARY- 
LAYER INTERACTION 

• ADDITIONAL VALIDATION EXPERIMENTS 

- HIGH SPEED MIXING 

- TRANSITION DUCTS 

- VORTEX GENERATORS 

t 
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EXPERIMENTS IN 3-D FLUID MECHANICS 


HYPERMIXING 



TRAN, DUCT 


VORTEX GEN 



k Yort«x generator*. 


Plane 4 


Plane 3 
Plane 2 


Plane 1 


CROSSING SHOCK 
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SHOCKS 


QUANTITATIVE DATA 
QUALITATIVE DATA 


UNSTEADY SHOCK 
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EXPERIMENTAL AND COMPUTATIONAL SURFACE AND FLOW-FIELD 
RESULTS FOR AN ALL-BODY HYPERSONIC AIRCRAFT 


WILLIAM K. LOCKMAN 
NASA Ames Research Center 
Moffett Field, CA 94035-1000 


SCOTT L. LAWRENCE 
NASA Ames Research Center 
Moffett Field, CA 94035-1000 

JOSEPH W. CLEARY 
Eloret Institute 
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OUTLINE 


OBJECTIVE 

EXPERIMENT 

• Model 

• Test Conditions 

• Measurements 

UPS CODE (Upwind PNS Solver) 

RESULTS (Experimental & Computational) 

• Flow Visualization 

• Surface Pressures 

• Surface Convective Heat Transfer 

• Pitot-Pressure Flow-Field Surveys 

CONCLUDING REMARKS 


U7 


ALL-BODY HYPERSONIC TEST PROGRAM 
FOR CFD CODE VALIDATION 


OBJECTIVE: Establish benchmark experimental data base for generic 

hypersonic vehicle shape for validation and/or calibration 
of advanced CFD computer codes 


MOTIVATION: Need for extensive hypersonic data to fully validate CFD 

codes to be used for NASP & other hypersonic vehicles 


APPROACH: Conduct comprehensive test program for generic all-body 

hypersonic aircraft model in Ames 3.5-ft Hypersonic Wind 
Tunnel to obtain pertinent surface and flow-field data over 
broad range of test conditions 


ALL-BODY HYPERSONIC AIRCRAFT EXPERIMENT 

Ames 3.5-ft HWT 


MODEL : 

AMES GENERIC ALL-BODY HYPERSONIC MODEL: 

• Delta Planform (A = 75°) 

• Forebody — Elliptic Cone (a/b = 4) 

• Afterbody — Elliptic Cross Sections with Sharp Trailing Edge 

• Sharp or Blunt Nose Tip 

• With or Without Control Surfaces (Tested without to date) 

• Canard 

• Combination Horizontal/Vertical Tails 


ALL-BODY HYPERSONIC AIRCRAFT MODEL 

W/O CONTROL SURFACES 


Elliptical 
Cross Sections 



Forebody - Elliptic Cone (a/b = 4) with Sharp or Blunt Nose Tip 
Afterbody - Elliptical Cross Sections with Sharp Trailing Edge 


ALL-BODY HYPERSONIC AIRCRAFT MODEL 
IN NASA/AMES 3.5 -FT HWT 

W/O CONTROL SURFACES; LENGTH = 3 FT 



ALL-BODY HYPERSONIC AIRCRAFT EXPERIMENT 

Ames 3.5-ft HWT 

TEST CONDITIONS : 

• M m = 5, 7, & 10 

• Re . = 1.5x10 s to25x10 6 (Laminar to Turbulent Flows) 

• a = 0° to 1 5° (Attached & Separated Flows) 

MEASUREMENTS : 

• Flow Visualization 

• Shadowgraphs 

• Surface Oil-Flow Patterns (Skin-Friction Lines) 

• Surface Pressures 

• Surface Convective Heat Transfer (Selected Areas) 

• Flow-Field Surveys 

• Pitot-Pressure Probes 

• Laser Doppler Velocimetry (To be done) 

— Mean Velocities 

— Turbulence Quantities 


UPS CODE 

(Upwind PNS Solver) 


CHARACTERISTICS OF ALGORITHM: 

• Second-order accurate and upwind in crossflow directions 

• First-order accurate in streamwise (marching) direction 

• Implicit 

• Finite Volume 

ADVANTAGES OF UPWIND SCHEMES: 

• Shock waves are captured sharply and without oscillation 

• User specification of smoothing parameters is not required 

PRESENT ASSUMPTIONS: 

• Laminar flow 

• Turbulent flow 

— Boundary-layer transition specified 
— Baldwin-Lomax turbulence model 

• Perfect gas (Used for this study) 

• Equilibrium or nonequilibrium air 
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a UPS CODE; TURBULENT 

TANGENT-CONE METHOD; INVISCID 

TANGENT-WEDGE METHOD; INVISCID 


j ' i 1 

0 5 10 15 20 

LOWER-SURFACE INCLINATION ANGLE, 0 , deg 

W . ' 


/ 35 " 






SURFACE OIL-FLOW PATTERN 
a = 15°; M = 7.4; Re . =15 x 10 6 

00 °° j L 

LEEWARD 



PARTICLE TRACES FOR LEEWARD SURFACE BY UPS CODE 
(Simulated Surface Oil-flow Patterns; x/L = 0.0 to 0.9) 

a = 15°; M = 7.4; Re . = 15 x 10 6 

OO OOjL 
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EFFECT OF ANGLE OF ATTACK ON 
CENTERLINE SURFACE PRESSURES 
M„ = 7.4; Re„ L = 15 x 106 


WINDWARD 







j 
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EFFECT OF ANGLE OF ATTACK ON 
CENTERLINE SURFACE PRESSURES 
= 7.4; Re TOjL = 15 x 106 

LEEWARD 



AXIAL STATION, x/L 





SPANWISE SURFACE-PRESSURE DISTRIBUTIONS FOR 

FOREBODY AND AFTERBODY 

a = 0°; IVL = 7.4; Re^L = 15 x 10 6 




SPANWISE SURFACE-PRESSURE DISTRIBUTIONS FOR 

FOREBODY AND AFTERBODY 

a = 15°; WL = 7.4; Re^L = 15 x 10 6 


WINDWARD 
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SPANWISE SURFACE-PRESSURE DISTRIBUTIONS FOR 

FOREBODY AND AFTERBODY 

a = 15°; M„ = 7.4; Re„, L = 15 x 106 

! ’ 

LEEWARD 






INVISCID METHOD 

TANGENT-CONE 

----- SWEPT-CYLINDER 








NORMALIZED HEAT-TRANSFER RATE, q/q 


EFFECT OF REYNOLDS NUMBER ON 
CENTERLINE HEAT TRANSFER 
- 7.4; H w /H t = 0.4 



AXIAL STATION, x/L AXIAL STATION, x/L 


q = Stagnation-point heat-transfer rate for reference sphere (R = 0.01 L) 
' s 





STANTON NUMBER, St 


CENTERLINE HEAT TRANSFER 
a = 0°; WL = 7.4; Re„ x = 15 x 106; H w /H t = 0.4 
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STANTON N 


CENTERLINE HEAT TRANSFER 
WL = 7.4; Re^ L = 15 x 106; H w /H 


St = q/((pu)JH t - H w )) 


I 


WINDWARD 


» — gan g □ □ a an 

LEEWARD 
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EXPERIMENT 

o WINDWARD 
□ LEEWARD 


UPS CODE 
— TURBULENT 


AXIAL STATION, x/L 
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CENTERLINE HEAT TRANSFER 
; IVL = 7.4; Re^ , = 15 x 106; H w /H 


St = q/((pu)JH t - HJ) 


WINDWARD 



LEEWARD 
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EXPERIMENT UPS CODE 

o WINDWARD TURBULENT 

o LEEWARD 


AXIAL STATION, x/L 






CENTERLINE HEAT TRANSFER 
a = 15°; = 7.4; Re^L = 15 x 106; H w /H t = 0.4 


St = q/((pu)JH t - H w )) 


WINDWARD 


® e s a g g fi 


LEEWARD 
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EXPERIMENT 
o WINDWARD 
□ LEEWARD 

© WINDWARD 


UPS CODE 
TURBULENT 


(Re . = 25 x 10°) 

°°jL 


AXIAL STATION, x/L 
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STANTON NUMBER, St 


EFFECT OF ANGLE OF ATTACK ON 
CENTERLINE HEAT TRANSFER 
NL = 7.4; Re„ L = 15 x 10«*; H w /H t = 0.4 





STANTON NUMBER, St 


SPANWISE HEAT-TRANSFER DISTRIBUTIONS FOR FOREBODY 
a = 0°; WL = 7.4; Re TO>L = 15 x 106; h w /h, = 0.4 
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STANTON NUMBER, St 


SPANWISE HEAT-TRANSFER DISTRIBUTIONS FOR FOREBODY 
a = 5°; IVL = 7.4; Re^L = 15 x 106; H w /H t = 0.4 


St = q/(( P u )JH t - H w )) 


2/3 L— 1 

L- 


x/L EXPERIMENT UPS CODE, TURBULENT 
0.5 o l 

0.6 O /' 


WINDWARD 


LEEWARD 


.4 .6 

SPANWISE STATION, y/y. 



STANTON NUMBER, St 


SPANWISE HEAT-TRANSFER DISTRIBUTIONS FOR FOREBODY 
a = 10°; WL = 7.4; Re M>L = 15 x 10 6 ; H w /H t = 0.4 


4 x 10 


St = q/((pu)JH t - H w )) 


2/3 L— I 
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x/L EXPERIMENT UPS CODE, TURBULENT 

0.5 o 

0.6 □ 


WINDWARD 


LEE WARD 

i 0 i • 

.4 .6 

SPANWISE STATION, y/y LE 
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STANTON NUMBER, St 


SPANWISE HEAT-TRANSFER DISTRIBUTIONS FOR FOREBODY 
a = 15°; = 7.4; Re„, |L = 15 x 1Q6; H w /H t = 0.4 



.4 .6 
SPANWISE STATION, y/y. 






COMPARISONS OF EXPERIMENTAL PITOT-PRESSURE DISTRIBUTIONS 

WITH INVISCID MODEL OF FLOW 
Afterbody Centerline; x/L = 0.8; M„, = 7.4; Re^ L = 15 x 106 










mm 






NORMALIZED DISTANCE FROM SURFACE, 


EXPERIMENTAL AND COMPUTATIONAL PITOT-PRESSURE 
DISTRIBUTIONS AT VARIOUS AFTERBODY STATIONS 

Centerline; a = 0°; M„ = 7.4; Re M) L = 15 x 10 6 


ENTIRE SHOCK LAYER 


EXPERIMENT; x/L UPS CODE 

o 0.7 TURBULENT 

o 0.8 
.15 - ° 0.9 

BOUNDARY-LAYER EDGE 

(Shadowgraph) 


PITOT-PRESSURE RATIO, p. „/(p. J 




NORMALIZED DISTANCE FROM SURFACE, z g /L 


EXPERIMENTAL AND COMPUT ATI^^ PITO^PRESSURE 
DISTRIBUTIONS AT VARIOUS AFTERBODY STATIONS 

Centerline: a = 0°; = 7.4; Re^ i = 15 x 10® 


NEAR WALL 




NORMALIZED DISTANCE FROM SURFACE, z /L 


EFFECT OF ANGLE OF ATTACK ON 
PITOT-PRESSURE DISTRIBUTIONS AT AFTERBODY STATION 
Centerline; x/L = 0.8; M«, = 7.4; Re„ jL = 15 x 10 6 


ENTIRE SHOCK LAYER 



0 .5 1.0 1.5 2.0 2.5 3.0 3.5 4.0 

PITOT-PRESSURE RATIO, p, _/(p. „) 

t,Z r t,2 oo 


/ £ 


NORMALIZED DISTANCE FROM SURFACE, zj L 


EFFECT OF ANGLE OF ATTACK ON 
PITOT-PRESSURE DISTRIBUTIONS AT AFTERBODY STATION 
Centerline; x/L = 0.8; NL = 7.4; Re M> L = 15 x 10 6 


NEAR WALL 



PITOT-PRESSURE RATIO, p t 2 /(p ( g)^ 




CONCLUDING REMARKS 

DESCRIBED EXPT. IN AMES 3.5-FT HWT WITH ALL-BODY MODEL 
OUTLINED UPS CODE (Upwind PNS Solver) 

PRESENTED EXPERIMENTAL & COMPUTATIONAL RESULTS: 

• Flow Visualization (Shadowgraphs & Oil-Flow Patterns) 

• Surface Pressure Distributions 

• Surface Convective Heat-Transfer Distributions 

• Afterbody Pitot-Pressure Surveys 

OBSERVATIONS MADE: 

• Significant changes from forebody (conical) to afterbody (nonconical) flows 

• Complex leeward flow at angle of attack with cross-flow separation and vortices 

• Generally good agreement between experimental and UPS code results for: 

• Shock-Wave Angles 

• Surface Pressures (some differences at higher angles of attack and near leading 
edge) 

• Surface Heat Transfer (some differences for afterbody and leeward flows) 

• Afterbody Pitot-Pressure Surveys (some differences for leeward flow at higher 
angles of attack and for inner region of viscous layer) 
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COMPARISON OF EXPERIMENT WITH CALCULATIONS USING 
CURVATURE-CORRECTED ZERO AND TWO EQUATION 
TURBULENCE MODELS FOR A TWO-DIMENSIONAL U-DUCT 
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Abstract 

In this paper experimental measurements are com- 
pared with Navier-Stokes calculations using seven dif- 
ferent turbulence models for the internal flow in a two- 
dimensional U-duct. The configuration is represen- 
tative of many internal flows of engineering interest 
that experience strong curvature, such as that in the 
turnaround duct of the Space Shuttle Main Engine 
powerhead. A previous paper showed that application 
of a simple mixing length model to predict this flow 
gives poor agreement with important features of the 
experiment. In an effort to improve the agreement, 
this paper tests several versions of the two-equation 
k - e turbulence model including the standard ver- 
sion, an extended version with a production range time 
scale, and a version that includes curvature time scales. 
Each are tested in their high and low Reynolds num- 
ber formulations. Calculations using these new mod- 
els and the original mixing length model are compared 
here with measurements of mean and turbulence ve- 
locities. static pressure and skin friction in the U-duct 
at two Reynolds numbers. The comparisons show that 
only the low Reynolds number version of the extended 
k-e model does a reasonable job of predicting the im- 
portant features of this flow at both Reynolds numbers 
tested. 

Nomenclature 

A + = 26 

AR = aspect ratio 

Cf = skin friction coeff. = r w /q re f 
C p = static pressure coeff. = (p - p T ef)/qref 
H = channel height 

k = turbulent kinetic energy 

ss §(< u' > 2 + < v' > 2 ) 

l = mixing length 

* Research Scientist. Member AIAA 
t Research Scientist 

t Supervisory Aerospace Engineer. Member AIAA 


M = Mach number 

p = static pressure 

p t = total pressure 

q = dynamic pressure 

R c = streamline radius of curvature 

Re = Reynolds number based on H and U re f 

r = radial dist. from center of curvature 

s = down. dist. from channel entrance on duct c/l 

U, V = longitudinal, vertical mean velocities 

u, v = longitudinal, vertical inst. velocities 

u',r' = long., vert. inst. turb. vel. fluctuations 

uV = inst. turbulent Reynolds stress 

u T = wall friction velocity = y/r w /p 

= dim. velocity = U/u T 
x, y = x, y coordinates 

yi = distance from wall i in mixing length equation 

y t + = dim. wall variable for wall i = yiU T jv 

5 = boundary layer thickness 

e = turbulent energy dissipation rate 

0 = angle into bend measured from bend entrance 

#c = von Karman’s constant = 0.4 

i/ = fluid kinematic viscosity 

i/ t =s turbulent or eddy kinematic viscosity 

p — fluid density 

r t = turbulent shear stress = -p < uv > 

t* = local wall shear stress or skin friction 

w = mean flow vorticity 

<> = RMS time average 

Subscripts 

i , o = inner, outer walls 

r, 6 = radial, tangential dir. in cyl. coord. 

ref = ref. conditions 

s, n = parallel, normal dir. in stream, coord. 
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I. Introduction 

Computational fluid dynamics (CFD) is coming to 
play an increasingly major role in the initial design 
or verification of the internal flow in rocket engine 
components. For example, a three-dimensional in- 
compressible Navier-Stokes code (INS3D) was recently 
used to guide a possible redesign of the hot gas mani- 
fold for the Space Shuttle Main Engine (SSME) pow- 
erhead. Changing from a three-circular-duct to a 
new two-elliptical-duct configuration greatly improved 
flow through the new powerhead as confirmed by 
experiment, 1 - 2 

In spite of the above success, before CFD can be ap- 
plied widely and with sufficient confidence for future 
rocket engine design, the codes must be calibrated and 
verified by comparing them against well-documented 
experimental data. Most engine components oper- 
ate at very high pressures and flow Reynolds num- 
bers life's), so the flows are fully turbulent. Thus the 
accuracy of the codes is critically dependent on the 
turbulence models that they contain. Unfortunately, 
most turbulence models have been developed and veri- 
fied only for flows with very mild extra strain rates. In 
contrast, the flows in rocket engine components may be 
subjected to very large extra strain rates or other per- 
turbations. It is clear that knowledge on the structure 
of internal shear layers subject to such effects is cru- 
cial for successful numerical simulation of these flows. 
The accuracy of any turbulence model proposed for 
internal flow calculations should first be evaluated by 
comparing calculations using it with well-documented 
internal flow experiments containing one or more of 
the expected elements. 

To date, few experimental studies of such internal 
flows have been reported. There are at least two rea- 
sons for this. First, CFD has only recently been ap- 
plied to flows in rocket engines, so there has been little 
need for code verification experiments. Second, de- 
tailed measurements in rocket engines is usually quite 
difficult because of limited probe and/or optical access. 
One of the more widely-studied extra strain effects on 
internal flows is that arising from streamline curvature. 
There have been many such studies on the effects of 
mild curvature, and Bradshaw 3 gives a comprehensive 
review’ of them. However, information on the effects 
of strong curvature and other large extra strain rates 
is very much needed. 

Four recent studies have begun to examine both ex- 
perimentally and theoretically some of the important 
effects present in strongly-curved internal flows. These 
studies investigate the internal flow in 180° turnaround 
or U-ducts (TAD’s), where the radius ratio of bend 
centerline radius to duct height is of order unity. Such 
a geometry closely simulates many of the important 
features of the flow in the TAD of the SSME pow- 
erhead, such as the presence of strong curvature and 
pressure gradients, unsteady separation and interact- 


ing shear layers on opposite walls. 

In one study, Sharma et al. 4 measured the flow us- 
ing hot wires in an axisymmetric TAD air tunnel with 
M = 0.1 and Re = 10 5 . (For reference, the flow in 
the SSME TAD is at M = 0.1 and Re = 10 7 .) Chang 
and Kwak 5 calculated this flow using INS3D with a 
simple mixing length turbulence model. Both exper- 
iment and theory indicated a small amount of sepa- 
ration on the inner (convex) wall near the end of the 
bend. Poor agreement was found for the velocities 
near the outer (concave) wall and for the turbulent 
shear stresses throughout the bend, however. 

In a second study. Sandborn 6 measured the flow us- 
ing a one-component laser Doppler velocimeter (LDV) 
in a two-dimensional (2D) TAD water tunnel with 
U Te j zz 2 m/s. and Re = 7 x 10 4 to 5 x 10 5 . The results 
showed no inner wall separation at the lowest Re and 
greater separation than did Sharma et al. 4 at the high- 
est Re. Chen and Sandborn 7 calculated this flow using 
a Navier-Stokes (NS) solver and the two-equation Jb — e 
turbulence model both in its standard high Re form 8 , 
and also with a curvature correction. 9 Few conclusions 
about the accuracy of the models could be drawn, how- 
ever, since the calculations proceeded only to the 180° 
bend location. To that point the flow field is mainly 
pressure-driven and turbulence models are of less im- 
portance. 

In a third study, Avva et al. 10 calculated the TAD 
flow of Sandborn 6 at a Re = 9 x 10 4 using the standard 
k - e model in both its high and low Re versions. At 
that Re, Sandborn’* data showed only a small amount 
of separation on the inner wall at the bend exit. The 
low Re k — e model predicted this, whereas the high 
Re model predicted no separation. In contrast, and 
to the puzzlement of the authors, the high Re model 
predicted the measured static pressure very well and 
the low Re one underpredicted it. 

Finally, in a fourth study, Monson et al. 11 measured 
the flow using a two-component LDV in a 2D TAD air 
tunnel with ikf = 0.1 and Re = 10 5 and 10 6 . The mea- 
surements showed a small amount of inner wall separa- 
tion at the lower Re, and a much larger amount at the 
higher Re. The flow was calculated using INS3D and 
a simple mixing length model. Poor agreement was 
found with the experiment in several respects. The 
static pressure drop in the bend was badly underpre- 
dicted, no separation was predicted at the higher Re, 
and the computed velocities near both walls were in er- 
ror through most of the bend and further downstream. 
It appeared that the turbulence model was not repro- 
ducing the experiment very well, but no conclusions 
as to the cause of the failure could be drawn from the 
comparisons that were made. 

Because of the generally poor agreement found by 
the above studies between theory and experiment for 
TAD-type flows, it is clear that improved turbulence 
models will be required to accurately compute them. 
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ther downstream are much too low. For the outer wall 
(Fig. 9b), the model indicates separation at the bend 
entrance, whereas none was measured. The model pre- 
dicts very low Cj levels throughout the bend (as do all 
the models). Downstream Cj is also underestimated. 
The net result of the low Cj for the m.l. model is 
a prediction of too high a static pressure on the in- 
ner wall downstream of the bend (Fig. 10a). Note 
that the outer- wall pressure in the bend (Fig. 10b) is 
well-predicted by this and all of the other models. At 
Re = 10 6 . the m.l. model underpredicts Cj on both 
walls (Fig.'s 9c. 9d) to an even greater extent than at 
the lower Re. Also, attached flow is indicated, whereas 
the experiment shows increased separation. The net 
result is. once again, too high astatic pressure estimate 
downstream of the bend (Fig. 10c). It thus appears 
as though this zero-equation model, and others that 
employ the van Driest inner-scale mixing length (Eqn. 
3). are not well-suited to calculating this TAD-type of 
flow. 

Turning to the Cf results for the two standard k — e 
models at Re = 10 s in Fig. 9a. one can see that the low 
Re version (std. CH) shows early separation and later 
reattachment as compared to the experiment. (Note 
that the u waviness" of C f predicted by the CH model 
and other ITW models was also seen by Patel et al. 17 
for b.l. flow in the presence of adverse pressure gra- 
dients.) A large overshoot in Cj then occurs. This is 
typical behavior of the standard Chien low Re model 
in reattachment regions. 30 The high Re model predicts 
separation length better, but undershoots measured 
Cj further downstream. For the outer wall (Fig. 9b), 
the CH model agrees quite well with measured down- 
stream values, whereas the high Re model overpredicts 
the experiment. The net result for static pressure for 
these models (Fig. 10a), is that the CH version pre- 
dicts too low a pressure downstream of the bend, and 
the WF version agrees with experiment. For Re = 10®, 
both standard models indicate almost no separation 
(Fig.9c), which is opposite the experimental Re trend. 
Downstream of the bend, the WF version predicts Cj 
well on the inner wall, and the ITW version agrees 
with experiment on the outer wall. The net result 
is that both versions predict inner- wall pressure (Fig. 
10c) slightly higher than the measurements. It may 
be concluded that, although the standard k — e mod- 
els predict the TAD flow better than the m.l. model, 
neither version consistently predicts measured Cj , C v 
and extent of separation at both Re tested. Recent 
calculations (not included in this paper) of this TAD 
flow using the w LBl n17 low Re model rather than the 
CH model in std. k — e show much improved Cj and 
C v predictions. The prediction of separation extent is 
not improved, however. 

The next models to be considered for their Cj and 
Cp behavior are the two versions of the curvature- 
corrected k - e model. Recall that the goal of these 


models is to decrease TKE and TSS on the convex wall 
to be more consistent with experimental observations. 
It is obvious from the Cj data in Fig/s 9a and 9c that 
the models grossly overcorrect for such effects. Much 
too large an extent of separation is predicted at both 
Re. Following reattachment, both models also suffer 
from the same problems as did the standard models 
(i.e., an overshoot of Cj for the ITW version and an 
undershoot for the WF version). The net result is a 
prediction of too much downstream pressure loss for 
both versions at both Re (Fig.'s 10a and 10c). These 
results indicate that curvature corrections developed 
for flows with small curvature don’t always apply to 
flows with large curvature like the present TAD-type 
of flow. Thus, the curvature k - e model of Park and 
Chung 25 cannot be recommended for this flow. 

Finally, consider the skin friction and pressure pre- 
dictions of the extended k-e models. At Re = 10 5 , the 
ITW version indicates slightly more separation than 
the experiment Fig. 9a), but predicts the downstream 
recovery of Cj very well. The WF version shows much 
larger separation and an undershoot of downstream 
Cj. For the outer wall (Fig. 9b), the ITW model 
once again predicts downstream Cj and the WF model 
shows high values. For C p (Fig. 10a), the ITW model 
predicts pressure exactly and the WF model shows 
too much pressure loss. For Re = 10®, the ITW model 
predicts extent of separation and downstream Cj al- 
most exactly (Fig. 9c). (This is the only model tested 
that doesn’t show decreased separation at this Re.) 
The WF version once again shows very slow reattach- 
ment and downstream Cj recovery. Both versions of 
this model show excellent agreement with measured 
pressure at this Re (Fig.lOd). The conclusion can be 
reached that the extended fc — e model in its low Re ver- 
sion is the only one tested in this study that predicts 
most of the important measured features of this TAD- 
type of flow at both Re. The high Re version using 
the wall function formulation of Ref. 20 does not do as 
well, however. Perhaps a more sophisticated wall func- 
tion treatment of the separated region which includes a 
wake-like parameter 30 in the w law-of- the- wall”, and/or 
better choice of near- wall grid spacing *n certain re- 
gions of the TAD flow, would provide improved results. 
Overall, the extended k - e model together with the 
“LBl” low Re model overcomes the two main short- 
comings ot the standard CH k - e model for this TAD 
flow: 1) It reduces TKE and TSS on the convex wall 
so that separation is predicted at both Re; and 2) It 
eliminates the large Cj overshoot downstream of reat- 
tachment. Even though it is not developed specifically 
for curved flows, the extended LB k - e model seems 
to capture most of the important features of this flow 
and yet retains complete generality. Other turbulence 
models that reduce the level of turbulent stresses near 
the convex wall in the bend could perhaps do equally 
well when combined with the “LB1” low Re model. 


V. Summary and Conclusions 

A joint experimental and computational study of 

the flow in a two-dimensional U-duct with very strong 

curvature has been carried out. The significant con- 
clusions of this study are as follows: 

1* The experiment shows: a) significant turbulence 
enhancement on the outer (concave) wall consis- 
tent with a previously-observed curvature instabil- 
ity mechanism; b) almost total destruction of tur- 
bulence on the inner (convex) wall consistent with 
previous studies of flows with much less curvature: 
c) separation on the inner wall at the bend exit 
that increases with Reynolds number; and d) ex- 
treme core turbulence downstream of the bend that 
creates linear “plug-flow” type velocity profiles. 

2. The mixing length model of Patankar et al. 22 , and 
the standard and curvature-corrected k - e mod- 
els used with the low Reynolds number model of 
Chien lfr . failed to predict this flow and the trends 
with Reynolds number. 

3. The extended k - * turbulence model of Chen and 
Kim 24 , combined with the low Reynolds number 
model of Lam and Bremhorst 19 , was the only turbu- 
lence model of those evaluated that predicted mea- 
sured extent of separation, skin friction and static 
pressure throughout all regions of the flow at both 
low and high Reynolds number. The good predic- 
tion of skin friction downstream of the bend on the 
concave wall side may mean that modeling the de- 
tails of that complex flow in the bend itself may not 
be required. 

4 Application of wall functions as formulated by 
Launder and Spalding 20 to the above extended 
model produced poor results. More sophisticated 
treatment of the separated region and better choice 
of near- wall grid spacing in some regions of the flow 
may be required to improve predictions using wall 
functions. 
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Langley Investigation 


Instrumentation Techniques used: 

• 3-component Laser Velocimeter 

• 2-component hot wire 

• Pitot pressure surveys 

• 4th LV component for measuring freestream velocity 
Measurements to be made: 

• Mean Velocity Pitot, LV 

• Fluctuating Velocities LV, Hot wire 

• Spectral Information LV, Hot wire 

Measurement locations: 

• Global surveys - 700 pt surveys on grid (0.5 inch spacing) 

• Detailed surveys at x/H = -3, -2, -1, ..., 8, 9 (0.06 inch spacing) 

• Spectral data at x/H = 0, 7 

111 — July 1990 
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LV System Description 


Three-component Laser Velocimeter 

• 5 watt Argon ion laser 

• color separation 
(514.5,496.5,476.5 nm 
wavelengths) 

• sample volume size ~ 100 
microns 

• 0.8 micron polystyrene latex seed 
particle 

• data acquired in coincidence 

• orthogonal transmit optics 

• receive optics backscatter at 45 
degrees 

Fourth-component Laser Velocimeter 
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fiber-optic link using 488 nm 
wavelength 
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Mean Flowfield Characteristics 


Gross-flow Velocity Vectors 



Streamlines 
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Shear Stress Profile at x/H = 3 
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Concluding Remarks 

An investigation of the flow over a rearward facing step is currently 
underway. 

• A facility was constructed to meet the requirements of the 
DNS 

• Three-component LV data has been acquired throughout the 
flow documenting 

• Mean Velocities 

• Full Reynolds stress tensor 

• Preliminary review of the data by CFD personnel has 
expanded the scope to include additional Reynolds numbers 
and ER 

• Next phase of the test should begin Jan. 91 
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Plane 4 


Plane 3 
Plane 2 


Plane. 1 


VORTEX GEN 
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NASP-HYPERMIXING CONCEPTS 

EXPERIMENTAL INVESTIGATION IN 1x1 SWT 
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NASP-HYPERMIXING CONCEPTS 

TRACE GAS MEASUREMENTS 
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MASS FRACTION— MACH 1.6 

MASS FRACTION— MACH 3.0 


CD-89-44435 
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Vortex array embedded in a boundary layer. Top 
picture is the data taken at Mach .2. Bottom picture 
shows a model of the .data composed of Oseen vortices 
with vortex core locations , circulations , and peak 
Yorticities corresponding to the data results. 
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UNSTEADY SHOCK/BOUNDARY-LAYER 


UNIVERSITY OF TEXAS AT AUSTIN 
MACH 5 (DOLLING) 

PRINCETON UNIVERSITY 
MACH 3 (SMITS) 

LeRC 1X1 SWT 

MACH 3 AND 5 (BARNHART, HINGST) 
FLIGHT TESTING ? » 




Dynamics of Shock Wave 
Turbulent Boundary Layer Interactions 
Induced by Cylinders and Blunt Fins 

Principal 

Investigator: D. S. Dolling, Aerospace Engineering 

and Engineering Mechanics, 

The University of Texas at Austin 

Graduate Research Asst.: Leon Brusniak 

Sponsor: NASA Lewis Research Center 

(Warren Hingst) 
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Ensemble-Averaged Wall Pressure Distribution on 
Centerline Upstream of Cylinder (D = 1/2") 



Open Symbols: Shock-forward, y= Q-2, 

Closed Symbols: Shock-back, y= o- 8 


A \ 3 






Flow Conditions: 

Ma e = 2.9 
T 0 = 265 °K 
U e = 575 m/s 
Re e /m = 6.5 x 10^ 

Blunt fin diameter: 0.75" 

Boundary layer thickness: 1.1" 

Visualization: 

* Nd:YAG laser with a fourth harmonic generator operating at 266 nm 
(far UV), 4 ns pulses at 10 Hz. 

* Images obtained with a an intensified, UV sensitive CID camera. 

* Flow is from right to left 


Rayleigh scattering 

* In air flows: small ice clusters (30 nm) can dominate the Rayleigh 
signal. Yields high intensity signals in regions of low temperature, and 
signal dropout in regions of high temperature. A strong shock will 
appear as the boundary between a bright region of low density fluid and 
a dark region of high density fluid (the reverse of the expected result). 

* In nitrogen flows: the Rayleigh signal is now dominated by 
scattering from nitrogen molecules. A shock will appear as the boundary 
between a darker region of low density fluid and a brighter region of 
high density fluid (as expected). 


Rayleigh Scattering: 


* scattering by particles much smaller than the wavelength of light. 

* The intensity of the signal is proportional to the density of the 
scattering centers - gives the possibility of quantitative density maps. 

* Can reveal instantaneous turbulence and shock structure in a plane. 

* Sensitivity proportional to the fourth power of frequency, so 
operation at short wavelengths very beneficial. 


Blunt Fin Interaction Static Pressures uWvsz/v 

Sverdrup 
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NASA LOW SPEED CENTRIFUGAL COMPRESSOR 


• OBJECTIVES 

• COMPRESSOR PROGRAM 

• DESCRIPTION OF FACILITY AND COMPRESSOR 

• INSTRUMENTATION 

• CFD EFFORT & IMPACT ON RESEARCH PROGRAM 

• CFD/EXPERIMENT PRELIMINARY COMPARISONS 


SUMMARY 
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RESEARCH OBJECTIVES 


Improve understanding of flow in centrifugal compressors 


Provide data for modeling of various flow phenomena 


• Acquire "benchmark” data for 3-D viscous code validation 
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LOW SPEED CENTRIFUGAL COMPRESSOR 

BENCHMARK EXPERIMENTS FOR COMPUTATIONAL CODES 
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COMPRESSOR ASSEMBLY 
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NASA LOW SPEED CENTRIFUGAL COMPRESSOR 

DESIGN CHARACTERISTICS 


AIR FLOW 

66 lb m/sec 

ROTATIVE SPEED 

1920 RPM 

INLET PRESSURE 

14.7 psia 

PRESSURE RATIO 

1.18 

OUTLET TEMP. 

86 F 

EXHAUST 

ALTITUDE 
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NASA LOW SPEED CENTRIFUGAL COMPRESSOR 

INSTRUMENTATION 

• 465 STATIC PRESSURES (208 ON ROTOR) 

• 83 TOTAL PRESSURES (75 ON ROTOR) 

• 51 TOTAL TEMPERATURES (27 ON ROTOR) 

• 8 FLOW ANGLE PROBES 

• 28 STRAIN GAUGES 

• 6 PROXIMITY PROBES 


LASER ANEMOMETRY 


National Aeronautics and 

Space Administration INTERNAL FLUID MECHANICS DIVISION 

Lewis Research Center 
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LOW SPEED COMPRESSOR RESEARCH PROGRAM 

CENTRIFUGAL COMPRESSOR - CFD APPLICATIONS 

OBJECTIVE: 

Perform CFD analysis in order to provide guidance for 
experimental measurements: 

- develop a "feel" for the flow physics 

- aid in planning location/extent of measurements 

- on-line assessment of CFD limitations 


APPROACH: 

3D steady Navier-Stokes analysis of impeller only 

- VPI&SU grant 

- in-house 

3D steady Average -Passage analysis of impeller+diffuser 


STATUS: 

VPI&SU impeller analysis complete 
In-house impeller analysis in progress 
Average -passage stage analysis in progress 

LSCRP-3 





3D NAVIER STOKES CODE 

ANALYSIS OF LARGE LOW-SPEED 
CENTRIFUGAL IMPELLER 


PRESSURE DRIVEN SECONDARY FLOWS 
NEAR BLADE SUCTION SURFACE 


Blade leading 
edge 



- Blade trailing 
edge 






Axial Distance 



3D NAVIER STOKES CODE 
ANALYSIS OF LARGE LOW-SPEED 
CENTRIFUGAL IMPELLER 

PRESSURE DRIVEN 
TIP CLEARANCE FLOW 

Clearance 
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LOW SPEED CENTRIFUGAL COMPRESSOR 
DIFFUSER SHROUD 
66 Ibm/sec 
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LOW SPEED CENTRIFUGAL COMPRESSOR 

3-D NS CODE PREDICTION V.S. EXPERIMENT 


3-D NAVIER STOKES ANALYSIS 
PARTICLE TRACES 


IMPELLER EXIT 

RADIAL VELOCITY DISTRIBUTION 







LOW SPEED CENTRIFUGAL COMPRESSOR 

DIFFUSER HUB 

66 Ibm/sec 



3 5" SL. 
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LeRC COMPRESSOR FLOW PHYSICS 

EXPERIMENT SCHEDULE 
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Low Speed Compressor 
Facility - IFMD 

Centrifugal 


full stage 


Multistage Axial 


4-stage axial 
compressor 


Small High Speed 
Compressor Facility - PSD 


centrifugal 


3-stage axial 

compressor 


compressor 


LSA-363 
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SUMMARY 

• LOW SPEED CENTRIFUGAL COMPRESSOR FACILITY 

V Improve understanding of complex flows 

V Data for flow physics modeling 

V "Benchmark” data for code assessment 

• CFD ANALYSIS 

V Develop feel for flow physics 

V Aid in planning location/extent of measurements 

V On-line assessment of CFD limitations 

• HARDWARE MODIFICATIONS 
^ 2% Tip clearance 

V Diffuser wedge plate 

• MILESTONES 

V Laser anemometer surveys LSCC, mid Sept. 

V High speed centrifugal, LFA surveys, 1990 
^ Low speed axial compressor, May 1991 

V Impeller/Diffuser LSCC, 1993 
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SECOND NASA CFD VALIDATION WORKSHOP 
NASA LEWIS RESEARCH CENTER 


PREPARED BY: LISA W. GRIFFIN , ROBERT WILLIAMS, AND PAUL MCCONNAUGHEY 

COMPUTATIONAL FLUID DYNAMICS BRANCH 

AEROPHYSICS DIVISION 

STRUCTURES AND DYNAMICS LABORATORY 


OVERVIEW 


DUCT AND PIPE BEND FLOW CODE VALIDATION 
90° DUCT BEND 

• OBJECTIVES 

• APPLICATION CONFIGURATION 
0 CODE METHODOLOGY 

% RESULTS 

• COMPUTER RESOURCE REQUIREMENTS 

• CONCLUSIONS 

TURBOMACHINERY CODE VALIDATION 
TURBINE CASCADES 

• OBJECTIVES 

• APPROACH 

• APPLICATIONS CONFIGURATIONS 

• CODE METHODOLOGY 

• RESULTS 

• COMPUTER RESOURCE REQUIREMENTS 

• CONCLUSIONS 
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DUCT AND PIPE BEND FLOW CODE VALIDATION 


CODE 

APPLICATION 

PRINCIPAL 

INVESTIGATORS 

INS3D 

* 90° DUCT BEND 
S-SHAPED DUCT BEND 
180° TURNAROUND DUCT 

PAUL MCCONNAUGHEY 
JONI CORNELISON 

INS3DLU 

* 90° DUCT BEND 
90° PIPE BEND 

ROBERT WILLIAMS 

FDNS3D 

* 90° DUCT BEND 
90° PIPE BEND 

ROBERT WILLIAMS 


OBJECTIVE 


ASSESS CODE CAPABILITY TO PREDICT COMPLEX SECONDARY 
FLOWS AND ASSOCIATED AXIAL FLOW VARIATION 


u o 


APPLICATION CONFIGURATION 

• NINETY DEGREE DUCT BEND WITH CONSTANT, SQUARE CROSS SECTION 
(TAYLOR, WHITELAW, AND YANNESKIS) 

• RADIUS RATIO - 2.3 

• REYNOLDS NUMBER - 790 

• DEAN NUMBER - 368 

• LASER-DOPPLER VELOCIMETER DATA 

• STREAMWISE AND GAPWISE VELOCITY COMPONENTS 

• 8 STREAMWISE STATIONS 

• 50 MEASUREMENTS PER STATION 
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CODE METHODOLOGY 


• INS3D 

• INCOMPRESSIBLE NAVIER-STOKES CODE 

• BEAM-WARMING /BRILEY MCDONALD ALGORITHM 

• ARTIFICIAL COMPRESSIBILITY 

• 2ND OR 4TH ORDER NUMERICAL DISSIPATION 
MODEL 

• INS3DLU 

• INCOMPRESSIBLE NAVIER-STOKES CODE 

• LOWER-UPPER SYMMETR1C-GAUSS-SIDEL IMPLICIT SCHEME 

• ARTIFICIAL COMPRESSIBILITY 

• 3RD ORDER NUMERICAL DISSIPATION MODEL 

• FDNS3D 

• UNSTEADY REYNOLDS-AVERAGED NAVIER-STOKES CODE 

• PRESSURE-BASED PREDICTOR-CORRECTOR PISO TYPE ALGORITHM 

• VARIABLE ORDER UPWIND DISSIPATION 

• K-b TURBULENCE MODEL FOR TURBULENT FLOWS 

• EQUILIBRIUM/FINITE RATE CHEMISTRY FOR REACTIVE 
FLOWS 




Streamwlse Velocity Contours at 


3 Plane 0.25 H Past the 90° Bend Exit. 



Inner Wall 










Gapwlse Velocity Contours at the Plane 0.25 H Past the 90° Bend Exit. 




EX! 


c 


e = 60° 

Cont. Int. = 0.045 



X = +0.25 H 
Cont. Int. = 0.03 


Predicted Radial Flow In 90° Bend Con 1 



6 = 77.5° 
Cont Int. s 0.025 



X =+2.5 H 
Cont Int. = 0.02 





X = +0.25 H 


X = 2.5 H 


Predicted Axial Flow in 90° Bend Compared with the Data of Taylor et al 
Contour Interval of Plots Is 0.08. 






GAP WISE DISTANCE FROM DUCT CENTERLINE 


STREAMWISE VELOCITY PROFILES - INS3DLU 




22,500 GRID POINTS 


49,000 GRID POINTS 


85,000 GRID POINTS 






0.25 H Past the 90° Bend Exit 








GAPW1SE DISTANCE PROM DUCT CENTERLINE 


0.4 


0.2 


0.0 


- 0.2 
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GAPWISE VELOCITY PROFILES - INS3DLU 



22,500 GRID POINTS 



49,000 GRID POINTS 85,000 GRID POINTS 


0.25 H Past the 90° Bend Exit 
























GAPVISE DISTANCE FROM DUCT CENTERLINE 



















COMPUTER RESOURCE REQUIREMENTS 

* ALL CALCULATIONS WERE RUN ON THE CRAY XMP-416 AT MSFC 


• INS3D 

• 17 WORDS/GRID POINT OF MEMORY 

• 25.0 MICRO SEC/GRID POINT/ITERATION OF CPU TIME 

• REQUIRED 3.0 MW OF MEMORY AND 1.2 CPU HOURS FOR THE CALCULATION WITH THE 
DENSEST GRID IN THIS STUDY BASED ON 1000 ITERATIONS 


• INS3DLU 

4 

• 74 WORDS/GRID POINT OF MEMORY 

• 1 1 .5 MICRO SEC/GRID POINT/ITERATION OF CPU TIME 

• REQUIRED 6.1 MW OF MEMORY AND .14 CPU HOURS FOR THE CALCULATION WITH THE 
DENSEST GRID IN THIS STUDY BASED ON 500 ITERATIONS 

• FDNS3D 

• 46 WORDS/GRID POINT OF MEMORY 

• 127.4 MICRO SEC/GRID POINT/ITERATION OF CPU TIME 

• REQURIED 3.7 MW OF MEMORY AND 3.0 CPU HOURS FOR THE CALCULATION WITH THE 
DENSEST GRID IN THIS STUDY BASED ON 1000 ITERATIONS 



SUMMARY/CONCLUSIONS 


• VALIDATION OF INS3D, INS3DLU, AND FDNS3D FOR LAMINAR FLOW THROUGH A 90° DUCT 
BEND HAS BEEN COMPLETED 

• AGREEMENT BETWEEN PREDICTED AND MEASURED AXIAL FLOW WAS GOOD FOR EACH 
CODE (PROVIDED THE GRIDS WERE DENSE ENOUGH TO ADEQUATELY RESOLVE ALL FLOW 
COMPONENTS). 

• FOR ADEQUATE FLOW RESOLUTION, INS3DLU AND FDNS3D REQUIRED GRIDS OF 85,000 GRID 
POINTS. INS3D REQUIRED 100,000 TO 175,000 GRID POINTS. 

• INS3DLU SHOWS A MARKED IMPROVEMENT OVER INS3D IN COMPUTATIONAL SPEED. FDNS3D 
A MORE 

GENERAL CODE, IS THE LEAST EFFICIENT IN TERMS OF COMPUTATIONAL SPEED OF THE 
THREE CODES 


CODE 

MTSBL/STAN5 

FDNS3D 

ROTOR1 


ROTOR3 


TURBOMACHINERY CODE VALIDATION 


APPLICATION 

UTRC LSRR FIRST STAGE 
SSME HPFTP (FPL) FIRST STAGE 
SSME HPFTP TTA FIRST STAGE 
(MSFC AND CALSPAN TESTS) 


SSME HPFTP TTA FIRST STAGE 
(MSFC TESTS) 


UTRC LSRR FIRST STAGE 
SSME HPFTP (FPL) FIRST STAGE 

SSME HPFTP TTA FIRST STAGE 
(MSFC TESTS) 


KOPPER’S CASCADE 
HODSON’S CASCADE 

* UTRC LSRR FIRST STATOR 

* LANGSTON’S CASCADE 
GGGT FIRST STATOR 

SSME HPFTP TTA FIRST STAGE 
(CALSPAN TESTS) 

UTRC LSRR FIRST ROTOR 
HEAT TRANSFER 


PRINCIPAL INVESTIGATORS 

HELEN MCCONNAUGHEY 
HELEN MCCONNAUGHEY 
JOE RUF 


HELEN MCCONNAUGHEY 


LISA GRIFFIN 
LISA GRIFFIN, 

HELEN MCCONNAUGHEY 
LISA GRIFFIN 


LISA GRIFFIN/PRATT & WHITNEY 
LISA GRIFFIN/PRATT & WHITNEY 
LISA GRIFFIN/PRATT & WHITNEY 
LISA GRIFFIN/PRATT & WHITNEY 
LISA GRIFFIN 
LISA GRIFFIN 

LISA GRIFFIN 


INTRODUCTION/OBJECTIVES 


• OVERALL OBJECTIVE - VERIFICATION OF THE UNSTEADY, 3D, VISCOUS FLOW 

PREDICTION CAPABILITIES OF THE ROTOR3 CODE 

• PHASE I OBJECTIVE - VERIFICATION OF THE PERFORMANCE AND HEAT LOAD 

PREDICTION CAPABILITIES OF ROTOR3 THROUGH 
COMPARISON WITH BENCHMARK DATA 

• PHASE II OBJECTIVE - VERIFICATION OF THE UNSTEADY FLOW PREDICTION 

CAPABILITIES OF ROTOR3 THROUGH COMPARISON WITH 
UNSTEADY FLOW DATA 

• PHASE III OBJECTIVE - PREDICTION OF THE FLOWFIELD IN A ROCKET ENGINE 

TURBINE (SSME OR ADVANCED TURBINE) WITH ROTOR3 
AND COMPARISON OF THE PREDICTION WITH 
EXPERIMENTAL DATA 
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APPROACH 


• MODIFICATIONS WERE MADE TO ROTOR3 TO 

• ENHANCE THE CODE'S CAPABILITIES 

• DECREASE THE AMOUNT OF RUN TIME NECESSARY FOR PHASE I CODE 
VERIFICATION 

• CODE MODIFICATIONS 

• CREATION OF CASCADE VERSION OF ROTOR3 WITH WHICH TO VERIFY THE 
CODE'S PREDICTIVE CAPABILITIES IN STEADY FLOWS FOR WHICH DETAILED 
DATA IS AVAILABLE 

• INCORPORATION OF INLET BOUNDARY LAYER SPECIFICATION 

• INCORPORATION OF HEAT TRANSFER PREDICTION CAPABILITY 

• INCORPORATION OF TRANSITION LOCATION SPECIFICATION 


APPLICATION CONFIGURATION 


• UNITED TECHNOLOGIES RESEARCH CENTER LARGE SCALE ROTATING RIG (LSRR) 
FIRST STATOR 

• 3D ANNULAR CONFIGURATION 

• COMPARISON BETWEEN PREDICTIONS AND DATA USED TO VERIFY HEAT TRANSFER, 
SECONDARY FLOW, AND PERFORMANCE PREDICTION CAPABILITY 

• LANGSTON’S CASCADE 

• 3D PLANAR CONFIGURATION 

• TWO DIFFERENT INLET BOUNDARY LAYER PROFILES 

• COMPARISONS BETWEEN PREDICTIONS AND DATA USED TO ASSESS CODE 
SENSITIVITY TO INLET BOUNDARY CONDITIONS 




SOLUTION METHODOLOGY 

• 3D UNSTEADY, THIN-LAYER, COMPRESSIBLE NAVIER-STOKES 

• FACTORED, ITERATIVE, IMPLICIT ALGORITHM, THIRD-ORDER ACCURATE UPWIND 
DIFFERENCING SCHEME 

• MODIFIED BALDWIN-LOMAX TURBULENCE MODEL, FULLY TURBULENT FLOW ASSUMED 

t 

• ADIABATIC SURFACES ASSUMED 

• UNIFORM INLET CONDITIONS 

• OVERLAID AND PATCHED O- AND H-TYPE GRIDS, ROTOR H-GRID SLIDES PAST STATOR 
H-GRID AND OUTER CASING 



TOTAL PRESSURE LOSS CONTOURS DOWNSTEAM 
OF THE LSRR FIRST STATOR 


Casino 



COARSE GRID/PRESCRIBED INLET BOUNDARY LAYER 


REFINED GRID/PR ESCRIBED INLET BOUNDARY LAYER 



GAP AVERAGED TOTAL PRESSURE LOSS 
DOWNSTREAM OF THE LSSR FIRST STATOR 



0.0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1.0 

Z/b 





Stanton Numbers At The LSRR First Stator Midspan 



- 1.0 - 0.8 - 0.6 - 0.4 - 0.2 0.0 0.2 0.4 0.6 0.8 1.0 


X/C 




5 - 8899 - 0-233 



PREDICTED STANTON NUMBER CONTOURS 
FOR THE LSRR FIRST STATOR 


Hub 



PRESSURE SURFACE 


Effects Of 
Secondary Flow 


*8. Hub 


Leading 

edge 



Casing 


SUCTION SURFACE 


Trailing 

edge 




L; 


5 - 11004 - 0-291 















TOTAL PRESSURE LOSS CONTOURS DOWNSTREAM OF LANGSTON’S 
1 CASCADE, THIN INLET BOUNDARY LAYER 


MIDSPAN 



ENDWALL 





i 


TOTAL PRESSURE LOSS THROUGH LANGSTON’S CASCADE. THIN INLET 

BOUNDARY LAYER 





TOTAL PRESSURE LOSS THROUGH LANGSTON’S CASCADE, NOMINAL INLET 

BOUNDARY LAYER 



RESOURCE REQUIREMENTS 

• CALCULATIONS RUN ON THE CRAY X-MP 416 AND CRAY X-MP-280 

• CALCULATION REQUIREMENTS DEPENDENT UPON INLET CONDITION AND GRID 
DENSITY 


• LSRR CALCULATIONS STARTED FROM FREESTREAM 

• LSRR CALCULATION WITH REFINED GRID AND SPECIFIED INLET BOUNDARY LAYER 
PROFILE REQUIRED 20 CPU HOURS AND 4.0 X 10 6 WORDS OF CORE MEMORY PLUS 
6.4 X 10 6 WORDS OF SSD 

• LANGSTON’S CASCADE, THIN INLET BOUNDARY LAYER AND FINE GRID CASE, 
REQUIRED 4.0 X 10 6 WORDS OF CORE MEMORY AND 3.7 X 10 6 WORDS OF SSD. A 
CONVERGED SOLUTION REQUIRED 16 CPU HOURS WHEN STARTED FROM 
FREESTREAM CONDITIONS AND 7 CPU HOURS WHEN STARTED FROM AN EULER 
SOLUTION. 

• LANGSTON CASCADE, NOMINAL INLET BOUNDARY LAYER AND FINE GRID CASE, 
REQUIRED 4.0X10 6 WORDS OF CORE MEMORY PLUS 5.8 X10 6 WORDS OF SSD. A 
CONVERGED SOLUTION REQUIRED 12 CPU HOURS WHEN STARTED FROM AN EULER 
SOLUTION. 


2 ? o 


DISCUSSION SUMMARY 


• ACCURACY OF TOTAL PRESSURE LOSS PREDICTION TIED TO GRID DENSITY AND 
INLET BOUNDARY CONDITION 

■ OVERPREDICTION OF MIDSPAN LOSS FOR EACH CASE. OVERPREDICTION 
INCREASES AS SECONDARY FLOW INCREASES (TURBULENT MIXING?) 

- UNDERPREDICTION OF CASING BOUNDARY LAYER FOR ANNULAR CASCADE 
(CURVATURE EFFECTS?) 

• PREDICTED MIDSPAN HEAT TRANSFER SHOWS EXCELLENT AGREEMENT WITH 
DATA. COMPUTED STANTON NUMBER CONTOURS DISPLAY FEATURES CONSISTENT 
WITH THE FLOW 

i 

• HIGH RESOLUTION SMOOTH RESULTS TIED TO VERY FINE SPATIAL RESOLUTION 
AND, CONSEQUENTLY, LONG RUN TIMES. IMPROVED INITIALIZATION, SUCH AS 
STARTING THE CALCULATION WITH AN EULER SOLUTION, CAN REDUCE RUN TIME 
CONSIDERABLY. 


WORKING GROUP REGISTRATION 
LOW SPEED 
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Name 

Affiliation 

Gordon Pickett 

Pratt & Whitney 

Robert Simoneau 

NASA Lewis 

Scott Kjelgaard 

Nasa Langley 

Hugh Shin 

GE Aircraft Engines 

Don Dietrich 

GE Aircraft Engines 

Wei Tseng 

Naval Air Development Center 

Michael Crawford 

U. of Texas, Austin 

Terry Simon 

U. of Minnesota 

Mounir Ibrahim 

Cleveland State U. 

Ted Reyhner 

Boeing Commercial Airplanes 

Ed Schairer 

NASA Headquarters 

Mark Potapczuk 

NASA Lewis 

Lisa Griffin 

NASA Marshall 

Fred Gessner 

U. of Washington 

Reda Mankbadi 

NASA Lewis 

Michael Hathaway 

NASA Lewi s 
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